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Abstract 



We study jet physics in the high energy regime of QCD. Based on the 
NLO BFKL equation, we construct a vertex for the production of a jet at 
central rapidity in fc^-factorization. A jet algorithm is introduced, and we 
take special care of the separation of multi-Regge and quasi-multi-Regge 
kinematics. The connection with the energy scale of the evolution is inves- 
tigated in detail. The result is discussed for two situations: scattering of 
highly virtual photons, which requires a symmetric energy scale to separate 
the impact factors from the gluon Green's function, and hadron-hadron col- 
lisions, where a non-symmetric scale choice is needed. For the second case 
we are able to define a NLO unintegrated gluon density, valid in the small-x 
regime, and give the evolution equation for this density as well. 

In the second part, we examine the angular decorrelation of Mueller- 
Navelet jets. Using an operator formalism in the space of anomalous dimen- 
sion and conformal spin, we implement the NLO BFKL Green's function 
to study the rapidity dependence of angular decorrelations. We incorporate 
the necessary summation of collinearly enhanced corrections beyond NLO 
accuracy. We compare our results with data from the Tevatron pp-collider 
and provide predictions for the LHC as well. We also extend our study to 
the angular decorrelation between a forward jets and the electron in deep 
inelastic ep scattering. The angular decorrelation has not been measured in 
DIS so far, but we give theoretical results for this observable which already 
implement the experimental cuts. 



Zusammenfassung 



Wir untersuchen Jet-Physik im Hochenergie-Regime der QCD. Basierend 
auf der NLO BFKL Gleichung konstruieren wie einen Vertex fur die Produk- 
tion eines Jets in zentraler Rapiditat in fcr-Faktorisierung. Ein Jetalgorith- 
mus wird eingefiihrt, und wir verwenden besondere Sorgfalt auf die Trennung 
von Multi-Regge und Quasi-Multi-Regge Kinematik. Die Verbindung zur 
Energieskala der Evolution wird detailiert untersucht. Das Ergebnis wird 
fur zwei Situationen diskutiert: der Streuung von hoch virtuellen Photonen, 
welche eine symmetrische Energieskala erfordert, um die Impaktfaktoren 
von der Gluon-Greenfunktion zu separieren, sowie Hadron-Hadron Kollisio- 
nen, bei denen eine asymmetrische Wahl der Skala erforderlich ist. Fiir den 
zweiten Fall sind wir in der Lage, eine NLO unintegrierte Gluonendichte zu 
definieren, die im Bereich kleiner x giiltig ist, und geben die Evolutionsglei- 
chung fiir diese Dichte an. 

Im zweiten Teil betrachten wir die Winkeldekorrelation von Mueller- 
Navelet Jets. Unter Verwendung eines Operatorformalismus im Raum von 
anomaler Dimension und konformem Spin implementieren wir die NLO 
BFKL Greenfunktion, um die Rapiditatsabhangigkeit der Winkelkorrelation 
zu studieren. Wir beriicksichtigen die notwendige Summation von kollinear 
verstarkten Korrekturen jenseits der NLO Genauigkeit. Wir vergleichen un- 
sere Ergebnisse mit Daten vom pp-Beschleuniger TEVATRON und stellen 
Vorhersagen fiir den LHC zur Verfiigung. Zudem dehnen wir unsere Un- 
tersuchung auf die Winkelkorrelation zwischen einem Vorwartsjet und dem 
Elektron in tiefinelastischer ep-Streuung aus. Die Winkelabhangigkeit in 
diesem Kontext wurde noch nicht gemessen wurde, aber wir geben theo- 
retische Resultate fiir diese Observable an, die bereits die experimentellen 
Schnitte beriicksichtigen. 
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Chapter 1 

Introduction 



Today, Quantum Chromodynamics (QCD) has established itself as the ad- 
equate theory to describe the strong force between fermionic quarks and 
bosonic gluons. One of its striking features that led directly to its discov- 
ery as a candidate theory of the strong interactions is asymptotic freedom. 
Asymptotic freedom implies that at short distances, i.e. whenever a hard 
momentum scale is present, quarks behave as almost free particles and, 
hence, we can apply perturbation theory. The other side of the coin is con- 
finement, i.e. the coupling increases with the distance, and at large distances 
only bound states of quarks and gluons exist which cannot be described by 
means of a perturbative expansion. This interplay of perturbative versus 
nonperturbative physics, arising from a quite frugal appearing Lagrangian, 
makes QCD a challenging theory with a very rich phenomenology. 

For the soft scale sector of QCD we have to rely on purely phenomeno- 
logical models, QCD sum rules or on lattice calculations which have made 
impressing progress over the last decade. However, lattice calculations can- 
not answer all our question on QCD, and due to its dependence on vast 
computer resources lattice calculations are still restricted in its possibilities. 
Even though every hard scale process is accompanied by soft scale physics in 
a real experiment, it is possible to disentangle these parts using factorization 
theorems which allow to describe the hard part within perturbative QCD 
by a fixed order expansion in the coupling a s . These theorems are strictly 
proven for a limited number of processes, but have shown its usability in a 
large number of processes. 

A manifestation of these theorems is the well known collinear factor- 
ization in the context of a hadron involved in a collision process. Due to 
this factorization, cross sections of hadronic interactions can be written in 
terms of a process-dependent hard matrix element convoluted with univer- 
sal parton density functions which are described by the Dokshitzer-Gribov- 
Lipatov-Altarelli-Parisi (DGLAP) equation [H [2j O H] . Because of the strong 
ordering of virtuality in the DGLAP evolution, the virtualities of the par- 
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tons entering the hard interaction could be neglected compared with the 
large scale of the probe, i.e. they are considered as being collinear to the 
proton. 

Whenever one considers kinematic regions characterized by two differ- 
ent large scales, this fixed order calculation is not sufficient since logarithms 
of these two scales appear at each order of the perturbation series. If the 
two scales are ordered, the logarithm of the ratio of the two scales becomes 
large and compensates the smallness of the coupling. Therefore, these log- 
arithms have to be resummed to all orders to justify a perturbative treat- 
ment. One famous example is the case of high energy scattering with fixed 
momentum transfer. If the center of mass energy s is much larger than the 
momentum transfer \t\ - the so-called Regge asymptotics of the process - 
the gluon exchange in the crossed channel dominates and logarithms of the 
type [a s ln(s/|i|)] n have to be resummed. This is realized by the leading log- 
arithmic (LL) Balitsky-Fadin-Kuraev-Lipatov (BFKL) El El E] equation 
for the gluon Green's function describing the momentum exchange in the 
t-channel. 

In the context of hadron collisions such a two-scale situation is easily 
realized if the scale of the hard process is large compared to Aqcd ; but 
nevertheless small compared to the center of mass energy s. We are then 
in the region of small longitudinal momentum fractions x of the partons 
entering the hard matrix element. Therefor, there are no longer grounds for 
neglecting their transverse momenta fcy. It is believed that in this case a 
better descriptions is given by the BFKL evolution equation. Here large log- 
arithms of the form [a s ln(l/x)] n are taken into account. Another evolution 
equation which resums these type of logarithms is the Ciafaloni-Catani- 
Fiorani-Marchesini (CCFM) evolution equation P QUI CEH E] • Just as for 
DGLAP, it is possible to factorize the cross section into a a convolution of 
process-dependent hard matrix elements with universal parton distributions. 
But as the virtualities and transverse momenta are no longer ordered (as 
it is the case in DGLAP evolution), the matrix elements have to be taken 
off-shell and the convolution made also over transverse momenta with the 
so-called unintegrated parton densities. This factorization scheme is called 
kj -factorization [131 E] or semihard approach pT5| [16] . 

The BFKL equation as it stands is aimed to describe inclusive quantities. 
If we consider for instance deep inelastic scattering of a virtual photon on 
a proton the BFKL equation predicts the observed steep rise of the proton 
structure function ^(a?) at small Bjorken-x. For this observable - measured 
at the electron proton collider HERA - logarithms of the type [a s ln(l/x)] n 
have to be resummed. However, for F2 the hard scale of the photon has 
to be connected with the soft scale of the proton such that it is not clear, 
whether F<i can be completely described by pure BFKL dynamics. More- 
over, the existing HERA data can be described by conventional DGLAP 
evolution as well. But also the region of large x provides the opportunity 



9 



to observe BFKL dynamics. If one considers the production of forward jets 
instead of the proton structure function ^2(^)5 large logarithms of the form 
[a s ln(xBj/xFj)] n appear, with xbj being the Bjorken-rr of the photon and 
xfj the longitudinal momentum fraction of the forward jet. In the case of 
two jets produced in the very forward direction at a proton proton collision, 
logarithms of the form [a s ln(s/Q 2 )] n appear, where Q is the scale of the 
forward jets. Mueller and Navelet [T7] proposed this dijet cross section as 
an ideal observable to be described by the BFKL equation. 

One also has to deal with BFKL dynamics if one aims to understand 
actually the physics behind jet production. This knowledge is an essential 
ingredient in phenomenological studies at present and future colliders. At 
high center of mass energies the theoretical study of multijet events becomes 
an increasingly important task. In the context of collinear factorization the 
calculation of multijet production is complicated because of the large number 
of contributing diagrams. There is, however, a region of phase space where 
it is indeed possible to describe the production of a large number of jets: 
the Regge asymptotics (small-x region) of scattering amplitudes. If the jets 
are well separated in rapidity, the according matrix element factorizes with 
effective vertices for the jet production connected by a chain of i-channel 
Reggeons. 

It turned out to be very fruitful to investigate more exclusive observables 
then the total cross section. Considering again the case of two forward jets 
produced in a hadron hadron collision, not only the cross section grows with 
increasing rapidity distance between both jets, but the additional emission 
of gluons between the two jets tagged in the Mueller-Navelet scenario lead 
to an angular decorrelation between those two jets. Since the angular decor- 
relation stems from the additional emissions between the colliding particles, 
this decorrelation should be also observable in other reactions which can be 
described by the BFKL equation. 

One of the corner stones of the BFKL approach is the Reggeization of 
the gluon, i.e. a process involving color octet exchange in the t-channel 
is described at lowest order by the exchange of one gluon, but asymptot- 
ically behaves as ~ (s/ty^K Furthermore, Reggeization implies that the 
trajectory j(t) = 1 + u>(t) passes through 1 (the spin of the gluon) at t = 
(the mass of the gluon). Although the trajectory itself is not infrared fi- 
nite, the divergences are canceled if real emissions, using gauge invariant 
Reggeon-Reggeon-gluon vertices, are taken into account. It is then possi- 
ble to describe scattering amplitudes with any number of well separated 
particles (jets) in the final state. The (a s lns) n resummation is known as 
leading-order (LO) approximation and provides a simple picture of the un- 
derlying physics. Nevertheless, it suffers from some drawbacks. One of them 
is the complete indetermination of the energy scale sq scaling the energy s 
in the resummed logarithms. Another handicap concerns the coupling a s 
which, at LO, is just a global parameter with its scale being not restricted 
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as well. These limitations can be removed if the accuracy in the calculation 
is increased, and next-to-leading (NLO) terms of the form a s (a s lns) n are 
taken into account [18} 119] . Diagrams of higher order then include those 
which contribute to the running of the coupling, and the correct matching 
of the different higher order contributions fixes the energy scale sq. 

While at LO the only emission vertex - the Reggeon-Reggeon-gluon ver- 
tex - can be identified with the production of one jet, at NLO also Reggeon- 
Reggeon-gluon-gluon and Reggeon-Reggeon-quark-antiquark vertices enter 
the game. The first part of this thesis is dedicated to the construction of an 
inclusive one-jet-production vertex from these different contributions. Due 
to these new emission vertices at NLO we have to introduce a jet definition 
discriminating between the production of one or two jets by two particles. 
It is not sufficient to simply start from the fully integrated emission vertex 
available in the literature [181 119] . Rather, we have to carefully separate 
all the different contributions in its unintegrated form first before we can 
combine them. By this procedure we will also be able to determine the right 
choice of energy scales relevant for the process. Particular attention is given 
to the separation of multi-Regge and quasi-multi-Regge kinematics. 

As it will turn out, the scale of the two projectiles in the scattering pro- 
cess has a large impact on the structure of the result. The jet vertex can 
not be constructed without properly defining the interface to the scattering 
objects. To show this, we will perform this study for two different cases: 
the jet production in the scattering of two photons with large and similar 
virtualities, and in hadron-hadron collisions. In the former case the cross 
section has a factorized form in terms of the photon impact factors and of 
the gluon Green's function which is valid in the Regge limit. In the latter 
case, since the momentum scale of the hadron is substantially lower than 
the typical k? entering the production vertex, the gluon Green's function for 
hadron-hadron collisions has a slightly different BFKL kernel which, in par- 
ticular, also incorporates some /^-evolution from the nonperturbative, and 
model dependent, proton impact factor to the perturbative jet production 
vertex. 

Our final expression for the cross section of the jet production in hadron- 
hadron scattering contains an unintegrated gluon density. This density de- 
pends on the longitudinal momentum fraction - as typical to the conven- 
tional collinear factorization - and on the transverse momentum kx- This 
scheme is known as /c-r-factorization and has been considered up to now only 
at LO. In fact our result, valid in the small- x limit, shows that it is possible 
to extend the /^-factorization to NLO. Nevertheless, we have to state that 
it is not obvious how to give a more general formulation of unintegrated 
parton densities at NLO for general x. 

One of the most famous testing ground for BFKL physics are the al- 
ready mentioned Mueller-Navelet jets |17j . The predicted powerlike rise of 
the cross section with increasing energy has been observed at the Tevatron 
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pp-collider [20], but the measurements revealed an even stronger rise then 
predicted by BFKL calculations. As another option, we direct our attention 
to a more exclusive observable within this process. In the second half of 
this thesis we study the azimuthal correlation between these jets. When 
we consider hadron-hadron scattering in the common parton model to de- 
scribe two jet production, we have the following picture in mind: a gluon 
or quark is struck out of each of the initial hadrons which then scatters off 
the other one and - after the process of hadronization - is observed as a jet 
of particles. In this back-to-back reaction we expect the azimuthal angles 
of the two jets always to be ir and hence completely correlated. But when 
we rise the rapidity difference between these jets, the phase space allows 
for more and more emissions leading to an angular decorrelation between 
the jets. In the academical limit of infinite rapidity, the angles should be 
completely uncorrelated. In the regime of large, but realizable rapidity dif- 
ferences the resummation of large logarithms calls for a description within 
the BFKL theory. Unfortunately, the leading logarithmic approximation 
[21\ [2"2~] overestimates this decorrelation by far. Improvements have been 
obtained by taking into account some corrections of higher order like the 
running of the coupling [23\ I24j . A full NLO description would incorporate 
the NLO Green's function - in fact not in its angular averaged form - and 
the NLO jet vertices for Mueller-Navelet jets [251126] . Due to the complexity 
of both pieces a fully analytical treatment seems cumbersome. Even on the 
numerical level investigations so far have been performed only for studies of 
the NLO kernel itself [22 [Ml [Ml ED] • 

We make an approximation to a full NLO calculation and study the an- 
gular correlation on the basis of the NLO kernel. To consider the jet vertices 
just at LO accuracy is justified by the fact that the angular correlation is 
mainly driven by the evolution kernel, especially if one considers the rapid- 
ity dependence of the correlation. Nevertheless, we will discuss the correct 
implementation of NLO jet vertices and its consequences. Using the frame- 
work of an operator formalism, developed in Refs. [HI [32l [33] , we study the 
important NLO feature, namely, that it determines the energy scale so and 
the running of the coupling. Furthermore, we carefully study the influence 
of different renormalization schemes and in this context also the need of a 
partial resummation of collinearly enhanced higher order terms. 

Since the angular correlation of Mueller-Navelet jets was measured at 
the Tevatron collider by the D0 collaboration [34] , we can compare our 
results with these data. We find that our NLO calculation including the 
resummation of collinearly enhanced terms improves the LO description 
significantly, but still misses the data. Furthermore, we formulate our results 
as a prediction for the forthcoming Large Hadron Collider. 

The studies of Mueller-Navelet jets always have been related to forward 
jet experiments in deep inelastic scattering (DIS). The analogies can be 
drawn for the angular decorrelation as well. Therefore, we transfer the 
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machinery developed in the framework of Mueller-Navelet jets to DIS. Some 
details change due to the different projectile, namely the electron, but the 
main features - as expected from theoretical side - appear to be the same 
as in the Mueller-Navelet case. Unfortunately, there are no experimental 
studies dealing with angular correlation in DIS, but we hope that our work 
draws the attention to this problem also from the experimental side. 

This thesis is organized as follows: In chapter [2] we introduce and ex- 
plain the BFKL equation at LO and NLO accuracy. We especially work 
out the different contributions to the NLO kernel and discuss its properties. 
Furthermore we elaborate on the resummation of terms beyond NLO accu- 
racy for general conformal spin which are necessary to stabilize the NLO 
kernel. Chapter [3] shows how to construct the inclusive jet production cross 
section from the NLO BFKL equation. We give an explicit formula for the 
jet vertex. Furthermore, we show how, in the context of hadron scattering, 
we are led to an NLO unintegrated gluon density, and give an evolution 
equation for this quantity. The discussion on angular decorrelation is con- 
tained in chapter 01 There, we derive analytic expressions for the angular 
correlation in proton-proton and electron-proton scattering. We perform a 
resummation of collinearly enhanced terms of higher order and compare our 
calculations with existing data. In chapter [5] we draw our conclusions. In 
appendix [X] we present an alternative subtraction term for the jet vertex, 
as derived in chapter [3] The appendix [B] comprises the technical details on 
the resummed kernel. 



Chapter 2 

BFKL equation 



This first chapter is dedicated to the derivation of the BFKL equation in 
leading and next-to-leading logarithmic approximation. Starting from its 
original derivation ^ EJ [7J [H] a number of different approaches to the BFKL 
equation has been presented. Our presentation is by no means extensive, 
but instead focuses on those points which will be most relevant for the rest 
of this study. In the first section we present the derivation of the leading 
order BFKL equation and in the second chapter we show its generalization 
to next-to-leading order. We take special care of the features, which are 
new compared to LO and, hence, specific for NLO. In the third section we 
discuss the properties of the NLO kernel of the BFKL equation. The fourth 
section contains an overview over the phenomenology of BFKL dynamics. 



Let us consider the case of the total cross section <tab in the scattering of 
two particles A and B. It is convenient to work with the Mellin transform 



acting on the center-of-mass energy s. The dependence on the scaling factor 
so belongs to the NLO approximation since the LO calculation is formally 
independent of sq. 

If we denote the matrix element for the transition A + B — ► A + B + 
n produced particles with momenta k{ (i = 1,... ,n) as -A^ +n , and the 
corresponding element of phase space as d&^ +n , we can write the total 
cross section as 



As we mentioned in the introduction we are interested in the Regge limit 
where s is asymptotically larger than any other scale in the scattering pro- 



2.1 BFKL equation at LO 




(2.1) 
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cess. In this region the scattering amplitudes are dominated by the produc- 
tion of partons widely separated in rapidity from each other. This particular 
configuration of phase space is known as multi-Regge kinematics (MRK). 
Particles produced in MRK are strongly ordered in rapidity, but there is 
no ordering of the transverse momenta which are only assumed not to be 
growing with energy. 

We fix our notation in Fig. I2.lt qi correspond to the momenta of those 
particles exchanged in the f-channel while the subenergies Si-n = (ki-i + 
ki) 2 are related to the rapidity difference between consecutive s-channel 
partons. Euclidean two-dimensional transverse momenta are denoted in 
bold. For future discussion we use the Sudakov decomposition ki = aipA + 
PiPB + ki± for the momenta of emitted particles. 




Figure 2.1: Notation for particle production in MRK. 



In MRK the center-of-mass energy for the incoming external particles 
can be expressed in terms of the internal subenergies as 



.1=1 



.1=1 



-1 
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q?q«+i II 

i=l 



Si-1, 



ktik? 



(2.3) 



where we have used the fact that in Regge kinematics s is much larger than 
—t and, therefore, Qo — @ n +l — 1, kg ~ and k^ +1 ~ q n+ i- To write 
down the measure of phase space we use dimensional regularization with 
D = 4 + 2e, i.e. 
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2.1. BFKL equation at LO 
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At tree level, the matrix element .A^B+n °^ Q2.2p can be written in MRK 
in the factorized form 



A- - 

AB+n 

2s 



n 1 



Qi+i) 



1 



a 2 



(2.5) 



with Tp being the couplings of the gluon to the external particles, and 
7(^,^+1) the gauge invariant effective (non-local) vertices of two t-channel 
gluons to a produced s-channel gluon. These vertices are called non-local 
since they encode the propagators of the connected t-channel gluons. 

A calculation of the loop corrections in leading logarithmic approxima- 
tion does not break the factorized form of the amplitude. Instead it leads 
to the Reggeization of the t-channel gluons, i.e. the gluon propagator of 
qi gets a multiplicative factor ( ) with uji = u>(qf) the gluon Regge 



trajectory depending on the momentum carried by the Reggeon. At LO the 
Reggeon scale sr is a free parameter. 

The statement that the gluon reggeizes means that the Reggeon carries 
the same quantum numbers as the gluon, that the so-called trajectory j(t) = 
1 + u(t) passes through 1 at t = 0, and that this Reggeon gives the leading 
contribution in each order of perturbation theory. The consistency of this 
picture demands that if we start from the assumption of Reggeization the 
calculation of all real and virtual corrections to the elastic scattering of 
two particles with color octet exchange in the ^-channel has to result in 
the exchange of simply one Reggeon. This consistency condition and other 
related ones are called bootstrap conditions and play an important role. They 
have been proven at LO [35] and NLO [36] as well. 

Hence we can write the matrix element -^B+n * n ^ s resummed form at 
leading logarithmic accuracy as 

— 1 fa. - 1 /a . , \ + 1 



AB+n 

2s 



1 

.i=i y * 



-i.i 



SR 
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in+l 



Sn,n+l 
SR 



T B , (2.6) 



Gathering all these elements together it is possible to write the Mellin 
transform of Eq. (|2.ip as the sum 



F(u J ,s ) = ^2^ n \u;,s ), (2.7) 

n=0 

with the contributions from the emission of n s-channel gluons being 
•^ {n) (^o) 
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The impact factors <3?p and the real emission kernel for Reggeon-Reggeon 
into a s-channel gluon KL r can be written in terms of the square of the 
vertices Tp and 7, respectively. The kernel fC r (q», qj+i) is defined such that 
it includes one gluon propagator on each side: (qfq 2 , A 1 . The integration 
over Sj-i^ in Eq. (|2.8|) takes place from a finite so to infinity. At LO terms of 
the form uj In k? or In sp can be neglected when the integrand is expanded 



in a s . Therefore, at this accuracy, Eq. (|2.8j) gives 



(27T) 
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where the poles in the complex w-plane correspond to Reggeon propaga- 
tors. This simple structure is a consequence of the linearity of the integral 
equation for the gluon Green's function. To see this connection explicitly 
we can introduce the gluon Green's function 



/u>(qa,q&) = ^ 



n=0 



" +1 d D -\, 



n 

i=i 



W - 2w; 



]f[/C r (qi,q i+ i) 



,i=l 



5 (2) (qa-qi)^ (2) (q6-q„+l) (2.10) 



being a solution of the BFKL equation 

oofMa^b) = <5 (2) (qa - q fe ) + J d D - 2 q/C(q a ,q)iUq,q fe ). (2.11) 

The integration kernel fC contains the real emission kernel /C r and the tra- 
jectory 

k->(q a ) which in this context often is called the virtual part of the 

kernel. 

/C(q a , q) = 2uj{c&) 8® (q a - q) + /C r (q a , q) (2.12) 

For the sake of clarity we have to say that Eq. (|2.1ip as derived here 
is only the forward BFKL equation. The name stems from the fact that 
the total cross section is related to the imaginary part of the elastic forward 
scattering amplitude by the optical theorem. The basis of the optical the- 
orem is unitarity of the S'-matrix which can be applied to the non forward 
scattering as well. The non-forward version of Eq. (|2,2p connects the imagi- 
nary part of the elastic scattering A + B — ► A' + B' with the sum of products 
of two different production amplitudes: 



5mi AB ^A'B' 



1 00 

-T 

n=0 



^AB+n^AB-^AB+n ("^A'B'^ABH 



V 
n/ 



(2.13) 



Following the same logic and steps leads to the BFKL equation with mo- 
mentum transfer. 
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Coming back to the case of the total cross section and summing the 
contributions of Eq. (|2.9|) we can write Eq. (|2.7[) as a convolution of this 
Green's function with the impact factors 

**(<«>. *>) = 75-4z2 // d °~l qa ^qa)/a,(q»,qt)<Mq&). (2-14) 

(2vr)^ ^ 

We will see below that Eq. (12, 9p holds very similarly at NLO. 

After this brief introduction to the structure of BFKL cross sections and 
its iterative expression we now turn to the NLO case. The factorization 
into impact factors and Green's function will remain, while the kernel and 
trajectory will be more complex than at LO. We discuss these points in the 
next section. 



2.2 Different contributions at NLO 



To discuss the various contributions to NLO BFKL cross sections we follow 
Ref. [37]. We comment in more detail those points which will turn out to 
be more relevant for our later discussion of inclusive jet production. Our 
starting point are Eqs. (12. ip to (12. 4p . which remain unchanged. Since at 
NLO the sr scale is no longer a free parameter, we should modify Eq. (|2.6p 
to read 



A 



AB+n 

2s 



l(Sfl;0,l) 
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.i=i Hi 
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7 (sfl;i - M ' SiiiM+l) (%%+i) 



■S_R;n,n+l 



-i( s fl;n,n+l) 



B 



(2.15) 



The propagation of a Reggeized gluon with momentum qi in MRK takes 
place between two emissions with momenta /cj_i and k{ (see Fig. l2.ip . There- 
fore, at NLO, the term sr, which scales the invariant energy Sj_i,i, does 
depend on these two consecutive emissions and, in general, will be written 
as SR-i-n. It is important to note that the production amplitudes should 
be independent of the energy scale chosen and, therefore, 



'i?;0,l 



for the particle-particle-Reggeon vertices and 



S R;n,n+l 



(2.16) 
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(2.17) 
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for the Reggeon-Reggeon-gluon production vertices. 

At NLO, besides the two-loop corrections to the gluon Regge trajectory 
|38] , there are four other contributions which affect the real emission vertex. 
The first one consists of virtual corrections to the one gluon production 
vertex |39} l4"0l I41j . The second stems from the fact that in a chain of 
emissions widely separated in rapidity two of them are allowed to be nearby 
in this variable which is known as quasi-multi-Regge kinematics (QMRK) 
[121 H3] . A third source is obtained by perturbatively expanding the Reggeon 
propagators in Eq. (j2. 15[) while keeping MRK and every vertex at LO. A 
final fourth contribution is that of the production of quark-antiquark pairs 
|44} [T3l [4"5l l4*6t ITf] . The common feature of all of these new NLO elements is 
that they generate an extra power in the coupling constant without building 
up a corresponding logarithm of energy so that a s (a s In s) n terms are taken 
into account. 

With the idea of introducing a jet definition later on, it is important to 
understand the properties of the production vertex which we will describe 
now in some detail. 

Let us u start with the virtual corrections to the single-gluon emission 
vertex. These are rather simple and correspond to Eq. (|2.9p with the in- 
sertion of a single kernel or impact factor with NLO virtual contributions 
(noted as (y)) while leaving the rest of the expression at Born level (written 
as (B)). More explicitly we have 
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(2.18) 



Now we turn to the discussion of how to define QMRK. For this purpose 
the introduction of an extra scale is mandatory in order to define a separation 
in rapidity space between different emissions. Following Ref. [37], we call 
this new scale s\. At LO MRK implies that all Sjj = {hi + kj) 2 are larger 
than s\. In rapidity space this means that their rapidity difference \yi — y~\ 

is larger than ln(,SAA/k 2 k 2 ). As stated earlier, in QMRK one single pair of 
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emissions is allowed to be close in rapidity. When any of these two emissions 
is one of the external particles A or B, it contributes as a real correction 
to the corresponding impact factor. If this is not the case it qualifies as a 
real correction to the kernel. This is summarized in the following expression 
where we denote real corrections to the impact factors by (r): 
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(2.19) 



The modifications due to QMRK belonging to the kernel or to the impact 

(r) 

factors are, respectively, /Cqmrk and <fr p , i.e. 



^QMRK(q 4 ,q J+ l =(^c - 1 / ds /„ N g 2 2 ' ( 2 - 20 ) 

' ( 27r ) q* q^+i 



*W(k) =^/^T / W^giy(£A^j) . (2 . 21) 



In both cases s denotes the invariant mass of the two emissions in QMRK. 
The Heaviside functions are used to separate the regions of phase space 
where the emissions are at a relative rapidity separation smaller than s\. It 
is within this region where the LO emission kernel is modified, grr^gg an d 
&pr-^pg are the total cross sections for scattering of two Reggeons into two 
gluons, and an external particle and a Reggeon into an external particle and 
a gluon, respectively. Irr and Ipr are the corresponding invariant fluxes, 
and N c is the number of colors. 

For those sectors remaining in the MRK we use a Heaviside function 
to ensure that Sj_i , > s\. In this way, MRK is clearly separated from 
QMRK. We then follow the same steps as at LO and use Eq. (12. 8p with the 
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modifications already introduced in Eq. (|2.15p . i.e. 
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After performing the integration over the Sj_i j variables the following in- 
teresting dependence on sa arises: 
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(2.23) 



It is now convenient to go back to Eq. (|2.ip and write the lower limit 
so of the Mellin transform as a generic product of two scales related to the 
external impact factors, i.e. sq = ^/sq-a sq-b- By expanding in a s the factors 
with powers in to and Wj it is then possible to identify the NLO terms: 
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To combine this expression with that of the QMRK contribution we should 
make a choice for sr. The most convenient one is SR-ij = ^Jsr-i sr-j, where 
for intermediate Reggeon propagation we use SR-i = k 2 , and for the connec- 
tion with the external particles sr-q = sq-a and s,R ; n+2 = sq ; b- We can then 
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write 
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This corresponds to the LO result for J r ( n + 1 ) plus additional terms where 
the uj — 2u>i factor cancels in such a way that they can be combined with the 
LO result of F {n \ 

The quark contribution can be included in a straightforward manner 
since between the quark-antiquark emissions there is no propagation of a 
Reggeized gluon. In this way one can simply write 
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The production kernel can be written as 
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(2.27) 



with ct rr ^qq being the total cross section for two Reggeons producing the 
quark-antiquark pair with an invariant mass s. 

The combination of all the NLO contributions together generates the 
following expression for the NLO cross section: 
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(2.28) 



where the NLO real emission kernel contains several terms: 

Mq», q»+i) = (4 B) + /cf LO) ) (q», q»+i) 
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(2.29) 
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with JCqq given by Eq. (|2.27p . The two gluon production kernel ICqg is 
the combination of /Cqmrk of Eq. (|2,2U[) and the MRK contribution in 
Eq. ([2725]) . It explicitly reads 



Irrvrr-^gg(s) AQa 
(2^q?qf +1 

2 



/C GG (qi,qm) = {N 2 c -l) f ds 
- f d D -\ 4 B )(q 4 ,q)4 B )(q,q m )ilnf 7 . (2.30) 

j 2 V (q» - q) (q«+i - q) / 

Below we will show that in the limit s\ — ► oo the second term of this expres- 
sion subtracts the logarithmic divergence of the first one. When computing 
the total cross section, it is natural to remove the dependence on the pa- 
rameter sa in this way. For our jet production cross section, however, we 
prefer to retain the dependence upon sa- 

For the impact factors a similar expression including virtual and MRK 
corrections as in Eq. (|2.2ip arises: 

<Mq W ) = 4 B) + + V^l / ^ /pfl ^t ) )? A "' ) 
J 2 v(qi - q) so;pJ 

From this expression it is now clear why to choose the factorized form so = 
y/so-A so;B- in this way each of the impact factors &a,b carry its own so-a,b 
term at NLO independently of the choice of scale for the other. 

To conclude this section, for the sake of clarity, the different contributions 
to the NLO BFKL kernel 
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Virtual 
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(b) 


QMRK 


n+1 


(c) 


MRK @ NLO 


n+1 


(d) 


Quark-antiquark pair 


n+1 


(e) 



are pictorially represented in Fig. 



As a final remark we would like to indicate that the divergences present 
in the gluon trajectories u){ (see Ref. [181 HH]) are au canceled inside the 
inclusive terms. We will see how the soft and collinear divergences of the 
production vertex are either canceled amongst its different components or 
are regularized by the jet definition. 

After having introduced the notation and highlighted the different con- 
stituents of a BFKL production kernel at NLO, we will discuss, in the coming 
section, the properties of this kernel in detail. 
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(a) (b) (c) (d) (e) 



Figure 2.2: Contributions to real emission kernel at LO (a) and NLO (b-e). 

2.3 Properties of the NLO BFKL kernel and re- 
summation of terms beyond NLO 

Shortly after the derivation of the NLO corrections to the BFKL kernel 
[TBI \19\ HU US] it was pointed out [501 EJJ that; these corrections are not only 
large and negative, but that they also lead to possible instabilities. After 
the suggestion that NNLO or even NNNLO terms might cure the problem 
[52| [53] , various attempts have been made to guess the higher order terms 
and to resum them [5J1 [55l [Ml E7J EHJ [59J [60] . Actually, these procedures 
have proven to be very successful. 

Although such a resummation does not influence the construction of an 
explicit jet vertex in the following chapter, it plays a vital role when we dis- 
cuss the impact of the gluon Green's function in chapter[U Therefore, we like 
to motivate the resummation by a study of the kernel itself. Furthermore, 
we will discuss its implementation and consequences. 

2.3.1 The LO kernel and its eigenfunctions 

The LO BFKL equation (|2.1ip is solved if one knows a complete set of 
eigenfunctions </> nj „(q). At LO these eigenfunctions are well known to be 

<^(q) = -±= (q 2 )^ e me , (2.32) 

where the transverse momentum is expressed in terms of its modulus y^q 2 " 
and its azimuthal angle 9. The eigenvalues associated with the eigenfunc- 
tions 0n,f(q) are uj(ti,v) = a s xo(\n\,l/2 + iv), where a s = a s N c /ir and 

X0 (n, 7 ) = 2^(1) -^( 7+ f)-^(l- 7+ ^). (2.33) 

The function ip is given as ip(x) = T'(x)/T(x), where T is the Euler gamma 
function. 

It is natural to change from the presentation in terms of transverse mo- 
menta to those of conformal spin n and anomalous dimension 7 = 1/2 + iv. 
We use this change to introduce an operator representation used in Ref. [31] 
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and extended in Ref. [32|, 133]. The transverse momentum representation is 
defined by 

qilqi) =qi|qi) (qi|q 2 ) =s {2) (qi - q 2 ) , (2.34) 

the kernel of the operator fC is 

£(qi,q 2 ) = (qi|/C|q 2 ), (2.35) 
and the BFKL equation together with its formal solution simply reads 

uL=L + Kf u = — (2-36) 

LO — K. 

As a second basis we will use n and v. The eigenfunctions given in 
Eq. (|2.32p relate both bases: 

( qi |^,n) = -l=(q?)^-5 e ^, (2.37) 

where 6\ is the azimuthal angle of qi . The normalization of this new basis 
follows directly as 

<n»,n> = J d^^^-^e^-^^S^-u^, (2.38) 
and the LO eigenvalue equation simply reads 

K.\v, n) = uj (n, v) \u, n) = a s xo ( \n\, - + if] \v,n). (2.39) 



Often the functions cj(n, v) and x( n )l/2 + iv) are called kernel as well. 
Knowing the action of the kernel (see Eq. (12.390 ). we can use the formal 
solution of the BFKL equation (|2.36p to write 

(n,v\LW,n) = z tttT r^8(v - v')6 nn >. (2.40) 

w - a s xo [\n\,2 

Using our operator formalism, it is straight forward to express the Green's 
function in physical terms of energy and transverse momenta: 

/ f qi,q2,— 
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For the last step we used Cauchy's theorem to perform the oj integration. 
Let us, now, focus on the conformal spin n = which is the dominant 
contribution. We will come back to higher conformal spins later. The 7 
dependence of the LO kernel is shown in Fig. 12. 31 It has a saddle point at 
7 = 1/2, and along the contour 7 = 1/2 + iv this saddle point appears as 
a maximum at v = 0. Since the kernel is exponentiated this pronounced 
maximum allows for a solution by a saddle point approximation, where one 
expands the kernel in v and keeps only the first two terms. 





Figure 2.3: LO kernel for conformal spin n = at a scale /i = 30GeV. 



2.3.2 The NLO kernel 

The action of the NLO kernel on the eigenfunctions of the LO one has been 
calculated in Ref. [61] and reads 



d 2 q 2 /C(qi,q 2 ) 



,in'(0 2 -0O 



a s (q?) 



X0 ,|n'|,i + ii/)+^)*(|n'|,i + «/ 



(2.42) 
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where 



<5(n,7) =4xi (n, 7) + 
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(2.43) 
(2.44) 



This means that the NLO kernel is diagonal with respect to the LO eigen- 
functions up to terms associated with the running of the coupling. While 
Xi( n 7l) is symmetric under the exchange 7 *-* 1 — 7, the remainder of 
S(n, 7) is antisymmetric. One could remove this part by modifying the LO 
eigenfunctions in the following way 
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Si 
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(2.45) 



Unfortunately, this set of functions no longer forms an orthonormal basis. 
Therefore, we stay with the LO eigenfunctions as they are. The antisymmet- 
ric term will vanish for some other reason in our study of angular decorre- 
lation, as we will show in Eq. (|4.24|) . The function xi f° r general conformal 
spin n is given by 
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(2.46) 



with the constant S = (4 - vr 2 + 5/3 /iV c )/12. CO) = EfcLi fc_n is the 
Riemann zeta function while the function 6 reads 
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Figure 2.4: NLO kernel for conformal spin n = at a scale [i = 30GeV. 

Let us have a closer look at the conformal spin n = component of the 
NLO kernel, depicted in Fig. 12.41 and compare it with the LO one, depicted 
in Fig. E3J 

The LO kernel has a single saddle point at 7 = 1/2. This property is 
changed in the NLO case if a s > .05. If we examine the contour 7 = 1/2 + iv 
for real v, we find, instead of a maximum at v = 0, two maxima at v ^ 
accompanied by a local minimum at v = 0. An unwelcome consequence of 
the saddle points having complex values is an oscillatory behavior in the 
ratio of the transverse momenta, which in the end can lead to negative cross 
sections. To understand the origin of this strange behavior, it is fruitful 
to investigate the pole structure of the kernel. The LO kernel (given in 
Eq. (|2.33[) ) has a simple pole at 7 = 0: 



In contrast, the NLO kernel (given in Eq. (12.460 ) has a more complicated 
pole structure around 7 = 0: 



The cubic pole does compensate for the equivalent terms appearing when 




(2.49) 




(2.50) 
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the symmetric energy scale is changed to an asymmetric one. The quadratic 
pole can be understood as a sequence of a small- x branching and a not- 
small- a; branching, since the explicit form of a turns out to equal 



A gg (0) + ^A qg (0) - (2.51) 



where A gg (co) {A qg (u))) is the Mellin transform of the splitting function 
Pgg(z) (P qg (z)). The /3q term reveals that also the running of the coupling 
generates a quadratic pole at 7 = 0. 

These higher-poles-terms are responsible for the problematic behavior of 
the kernel. The asymmetric scale change and the not-small- x branching are 
related to the collinear limit. Imposing that the collinear limit of the kernel 
does not conflict with the DGLAP equations, the problem can be cured 
\54\ [55] l56l [62] by a kernel which resums sub-leading collinear corrections. 
The focus, in these works, has been on conformal spin m = 0. Since in our 
study of decorrelation in chapter H] we are interested in angular dependences, 
we extend the solution proposed in Ref. [53] to general conformal spin. 

It can be shown \bA\ [5"5"] [5"6"l [6"2~] that this behavior is an artifact of the 
truncation of the perturbative expansion. If one starts from a proper form 
of the kernel in the DGLAP limit with an asymmetric energy scale, then the 
change to a symmetric energy scale for the BFKL regime can be obtained 
by a w-shift in 7. An expansion in u then gives the well known form of 
the NLO kernel. The resummation is based on the strategy to keep this u 
dependence in the kernel leading to a transcendental equation for the kernel. 

Technically the resummation is realized - for a given accuracy O (aj.) - 
by replacing divergences a"/ (7 + (n < I) with terms a™/(7 + y + f ) k 
without changing the correct expansion to O (a l s ). For a specific pole this 
change induces modifications of the order O (a"" 1 " 1 ) since u ~ a s \0 + O (o^) . 
Therefore, for all poles n < I we have to compensate this modification to 
keep the correct expansion to O (a l s ). Let us denote by x( m il) t ne BFKL 
kernel to all orders. Its expansion reads 

00 

xKt) = ^a"x n (m,7). (2.52) 

n=0 

So far the LO kernel xo an d the NLO kernel xi are known. The fixed order 
contribution as an expansion in 7 + ^ and 1 — 7 + y , respectively, reads 

Xn(m, 7 )=E 7 ih ^T + ° (7 + T ) (2.53a) 
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For the LO and NLO BFKL kernel we have d n k(m) = d n k(m) and 

d ol (m) =1, (2.54a) 
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2n f 



(2.54b) 

(2.54c) 
(2.54d) 



We now introduce a class of resummed kernels X (l)- Its perturbative 
expansion reproduces for the first N terms the expansion coefficients of the 
exact BFKL kernel. The differences to the exact kernel are of higher order. 



N 



x (A V 



7) = £*?Xn(m,7)+ £ «"^( m -7) ( 2 -55) 



n=0 



n=iV+l 



Although these resummed kernels will become uj dependent, we omit to 
write this dependence explicitly. Its expansion in 7 defines the coefficients 
^nk ( m ) m ^ ne following way 



00 d {N) (m\ 

k=o w + T ) 



( "V,t)=E it +^ + o(i- 7+ 



(2.56) 
(2.57) 



It is now necessary to define a set of functions £^(7) which are regular 
for 7 > 1/2 and behave at 7 = like l/7 fc . The concrete form of Dk(j) 
is of higher order accuracy. Developing the ideas presented in Ref. |54j . we 
construct from these functions a kernel X^°\l) which has its pole shifted by 
w/2: 



x (°)(m, 7 ) = xoKt) + do,iM [A + y + |) _jDi ( 7 + y 



+ do,i(m) 



fl 1-7 + t + t: -A 1 - 7 + 



7?) 



(2.58) 
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The master formula for N > then is a recursion formula and a generaliza- 
tion of Eq. (3.2) in Ref. 



X (JV) (m,7) = x {N 1} 0>7) + af (xN(m,-y) - xff ^(m,7)j 



fc=l 



dN,k(m) - d$ k (m) 



+ ( d Nik (m) -d)^, 1 ' (m) 



^(7+2 + 2)--°* ( 7+ y 



^(l-7+y+2)-^(l-7+2 



(2.59) 



The terms 



-AT (N-l) 



offdX?, , and —a^dX^, ^ are included to avoid 



double counting. The second and third line of Eq. (|2.59|) shifts the poles 
of xn by uj/2. The solution of the transcendental equation a s x^ (rn, 7) = ui 
for cj is the resummed kernel u(m, 7), which replaces a s Xo( m i ^)+a 2 s xi{rn, 7) 
Even though this procedure to modify the kernel is already quite con- 
crete, we want to keep the possibility to vary this procedure within its lim- 
itations. Consequently, we study different choices for ^(7) and do this in 
the same line as Ref. |54| . The explicit prescription of four different schemes 
can be found in appendix [Bj There we will discuss the different schemes 
which overall are very alike. Some minor differences let us favor scheme 3 
and we show the effect of the resummation, with scheme 3 as being represen- 
tative, in Fig. 12.51 The change to the pure NLO kernel is obvious: the poles 
at 7 = and 7 = 1 disappear and the shape along the contour 7 = 1/2 + iv 
has one single maximum at u = 0. 
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Figure 2.5: Resummed kernel (scheme 3) for conformal spin n 
fi = 30GeV. 



= at a scale 
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2.4 Phenomenology 

So far, we have introduced the BFKL equation at leading and next-to-leading 
order. Furthermore, we have discussed properties of the BFKL kernel and 
the need for a resummation of terms beyond NLO accuracy. Before we turn 
to consider specific observables in the following chapters, we would like to 
give some kind of overview over the phenomenology of BFKL dynamics. 

If one is willing to accept non perturbative input, a broad field of pro- 
cesses may be describable by the BFKL equation. We already mentioned 
the proton structure function F^, and diffractive scattering processes fall in 
this category as well. To study pure BFKL dynamics, it is necessary that 
both scattering partners provide a similar and hard scale. In the following 
we describe the phenomenology of processes that fulfill this necessity. 

One of the most famous predictions from the BFKL equation is the 
power-like rise s w of the total cross section with increasing energy. Due 
to lo being larger than 0, this behavior at some point violates unitarity, 
e.g. for hadronic processes the Froissart-Martin bound [63|. [M] restricts the 
growth of the total cross section like <r tot < Aln 2 s, where A is some con- 
stant. Therefore, for extremely high energies the descriptiveness by BFKL 
dynamics ends and different or modified models enter the game. 

But up to this regime, the power like rise should be observable, and the 
scattering of highly virtual photons has been proposed in Refs. [651 E3 EZ] as 
an appropriate environment to study this phenomenon. The virtual photons 
interact strongly when fluctuating into a quark-antiquark pair. The virtual- 
ities of the photons can be tuned to be of equal, or at least similar, scale such 
that DGLAP like evolution in the transverse momenta is suppressed offering 
a clean environment described by the BFKL equation. Therefore, it is often 
called the "gold-plated" process to test BFKL dynamics and is an object 
of intense studies from theoretical [681 E3 E3 [71] and experimental side 
[72J [73]. The large logarithm in this context is Y ~ In s 7 * 7 * / \J Q\Q\, where 
Q 2 and Q 2 are the photon virtualities. The data clearly show a rise with Y 
as expected by BFKL. Nevertheless, this rise is less steep than predicted by 
a LO BFKL calculation. 

Another candidate, summoned as a witness, are the already mentioned 
Mueller-Navelet jets. To be more explicit, we think of hadron-hadron colli- 
sions and of two jets emitted in the very forward/ backward region with a 
similar semihard scale p 2 : s S> p 2 ^> Aq CD . For large rapidity differences 
logarithms of the form [a s ln(s/p 2 )] n have to be resummed, which can be 
done in the framework of the BFKL equation. Mueller and Navelet proposed 
this process in Ref. [T7j as ideal to apply the BFKL equation and predicted 
a power like rise of the cross section. However, to see this growth of the 
cross section directly as a manifestation of Reggeization is hardly possible 
because it is drastically damped by the behavior of the parton distribution 
functions (PDFs) for x — > 1. One way out is to fix the PDFs and to vary the 
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center of mass energy of the hadron collider itself, and thereby vary the ra- 
pidity difference Ar] between the both jets. BFKL theory predicts a rapidity 
dependence of the cross section like a ~ exp((a — l)Ar/) jyj Ar], with a the 
so-called intercept . Usually, the collider energy is fixed and not tunable, 
but the D0 collaboration analyzed data taken at the Tevatron pp-collider 
from two periods of measurement. From these two points at = 1800GeV 
and 630GeV, they extracted an intercept of 1.65 ± .07 [20j . This is an even 
stronger growth of the cross section than predicted by a LO BFKL calcula- 
tion which for the kinematics of the D0 experiment yields a value of 1.45. 
It has been argued [73] that the exact experimental and theoretical defini- 
tions of the cross sections disagreed making an interpretation of the results 
difficult, and the fact that the experimental determination of the intercept 
is based on just two data points leaves room for more possible explanations. 

The measurement of forward jets at an electron-proton collider like 
HERA acts a combination of both processes and has been an object of 
intense studies as well. Compared to the proton-proton collider the depen- 
dence on the parton distribution functions can be singled out more easily. 
The large logarithm here is Y ~ Inxpj/xBj, where xfj is the longitudi- 
nal momentum fraction of the forward jet, and the variable XBj = Q 2 /s ep 
is given by the photon virtuality. By varying the virtuality of the photon 
one can tune XBj and thus tune Y without touching xpj- Measurements 
taken by the HI and the ZEUS collaboration [751 ITS'] have successfully been 
compared to BFKL calculations [TTJ EH] , while fixed order calculations un- 
derestimate the cross sections. Also the intercept has been extracted from 
these experiments |79] and is in better agreement with the theoretical pre- 
dictions then in the case of the Mueller-Navelet jets. Recently, new data on 
forward jet production at HERA have been published [801 El] ■ With the 
higher integrated luminosity, it was possible to present the triple differential 
cross section da/(dx^ dQ 2 dp^- et ), and first promising comparisons to NLO 
BFKL calculations have been presented [82] even though they rely on the 
saddle point approximation which for the NLO BFKL case is questionable 
as we will discuss in section 14.1.21 in detail. 

A very interesting experiment at an electron proton collider, as well, is 
the production of a vector meson - like p, <j) or J/tp - with a large momentum 
transfer —t » Aqcd- Considering the case of a large rapidity gap between 
the proton remanent and the vector meson, the assumption of just a two 
gluon exchange fails. This process was proposed in Ref. [83] as an observable 
to be described by the BFKL equation. In deed, further theoretical studies 
|84l l85j 186] are in very good agreement with the data from experiment [87] . 

For all these processes a complete NLO calculation is still outstanding. 
For the processes involving the virtual photon the reason is that the cal- 
culation of the NLO impact factor for the virtual photon is still work in 
progress (EH W\ EDJ EU E21 [93] . The NLO Mueller-Navelet jet vertices are 
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in principle available |25|, 126] . but are not yet cast in a form ready to be 
combined with the NLO Green's function. 

Instead, the impact factor for a virtual photon going into a light vector 
meson was calculated at NLO [M] up to terms suppressed by the virtuality 
of the photon. Quite recently, it has been used to calculate the electropro- 
duction of two light vector mesons in NLO accuracy within the framework 
of BFKL theory (31]. 

We will comment on the status for the angular observables in chapter [H 
where we present new calculations for these observables. Meanwhile we de- 
scribe - in the following chapter - how to calculate the inclusive production 
of jets in two different environments. The first one is the case of the interac- 
tion between two small and perturbative objects, namely highly virtual pho- 
tons, and the second will be the collision of two large and non-perturbative 
external particles such as the ones taking place at hadron-hadron colliders. 



Chapter 3 

Jet production 



In this chapter we study the inclusive production of a jet at central rapidity. 
We start with a LO description. With that first section we settle the no- 
tations and the framework, which we then - in the second section - extend 
to NLO. In doing so we show the necessary modifications if the scattering 
objects do not provide a hard scale. In the third section we discuss the 
separation between multi-Regge kinematics (MRK) and quasi-multi-Regge 
kinematics (QMRK) as well as the cancellation of divergences. 

3.1 Inclusive jet production at LO 

As MRK relies on the transverse scales of the emissions and internal lines 
being of the same order, it is natural to think that processes characterized 
by two large and similar transverse momenta are the ideal environment 
for BFKL dynamics to show up. Moreover, as the resummation is based on 
perturbative degrees of freedom, these large scales associated to the external 
particles should favor the accuracy of the predictions. An ideal scenario is 
the interaction between two photons with large virtualities Q\ 2 in the Regge 
limit s S> \t\ ~ Q\ ~ Q\. The total cross section for this process has been 
investigated in a large number of publications in recent years. Here we are 
interested in the inclusive production of a single jet in the central region 
of rapidity in this process. We will consider the case where the transverse 
momentum of the jet is of the same order as the virtualities of the photons. 

As a starting point we review single jet production at LO accuracy. As 
usual the total cross section can be written as a convolution of the photon 
impact factors with the gluon Green's function, i.e. 

. r d 2 K r d 2 k b r s+i °° du f s y 

(3-1) 

A common choice for the energy scale is sq = |k a | |k&| which naturally in- 
troduces the rapidities y^ and y^ of the emitted particles with momenta 
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and p g since 



v(ji&-Vb). (3.2) 

s o,' 



Let us remark that a change in this scale can be treated as a redefinition of 
the impact factors and, if so is chosen to depend only on k a or only on kj, 
the kernel as well. This treatment lies beyond LO and will be discussed in 
the next section. The gluon Green's function f u corresponds to the solution 
of the BFKL equation (pTTl) . 

For the inclusive production of a single jet we assign to it a rapidity 
yj and a transverse momentum kj, as shown in Fig. 13.11 In this way, if 

2 m fa- 



fe j = ckjpa + (3jPb + kj± the corresponding rapidity is yj = \ In tt- Using 



its on-shell condition we can write 




F^2 

PA + y^-e-y j pB + kj ± . (3.3) 






Figure 3.1: Total cross section and inclusive one jet production in the BFKL 
approach. 



It is possible to single out one gluon emission by extracting its emission 
probability from the BFKL kernel. The differential cross section in terms of 
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the jet variables can then be constructed in the following way: 



da 



d 2 kjdyj 



d 2 k a f d 2 k 



2vrk2 ./ 2vrkg 



$ 4 (k a )$ B (k 6 ) 

•<5+ioo 



x V(q Q ,qfe;kj,yj)^ 



' u - io ° dJ_ f.vu 
2iri \ s' 



/ w /(-q 6 ,-k 6 ) (3.4) 



with the LO emission vertex being 

V(q a , q fe ; kj, yj) = (q a , -q 6 ) 5^ (q a + q b - kj) 



(3.5) 



By selecting one emission to be exclusive we have factorized the gluon 
Green's function into two components. Each of them connects one of the 
external particles to the jet vertex. In the notation of Eq. (|3.4p the energies 
of these blocks are 



saj =(pA + qb) 2 , 



SB J =(pb + q a ) 2 - 



(3.6) 



In a symmetric situation, where both the jet and the impact factors provide 
a hard scale, a natural choice for the scales is similar to that in the total 
cross section 



so 



l k a| |k/| 



(3.7) 



These choices can now be related to the relative rapidity between the jet 
and the external particles. To set the ground for the NLO discussion of the 
next section we introduce an additional integration over the rapidity rj of 
the central system: 



da 



d 2 kjdyj 

d 2 K 



2nk 2 



d 2 q a / d 2 q b / drj 



S+iao j 
5—ioo ^ m 

5+ioo j, ,/ 



V(q a ,qfe,r/;kj,yj) 



dJ 
27d 



e u>> '(r\-y B ) f 



-q&, 



(3.f 



with the LO emission vertex being 

V(q a , qfe, m kj, yj) = IC { r B) (q a , -q b ) 5 (2) (q a + q& - kj) <5(r/ - yj). (3. 



Eqs. (|3.8p and (|3.9p will be the starting point for the NLO jet production 
in the symmetric configurations. 
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Let us now switch to the asymmetric case. In general we can write q a 
and q b as 

q a =a a PA + PaPB + Qa± Qb =a b p A + PbPB + Qb±- (3.10) 

The strong ordering in the rapidity of emissions translates into the conditions 
a a S> a b and /3& 3> /3 a . This, together with momentum conservation q a +q b = 
fej, leads us to aj = a a + a b « a a , (3j = f3 a + (3 b rs /3 6 and 

saj =A/s, s B j =ajs. (3.11) 

While the longitudinal momentum of q a (qb) is a linear combination of pa 
and pb we see that only its component along j?a(ps) matters. 

If the colliding external particles provide no perturbative scale as it is the 
case in hadron-hadron collisions, then the jet is the only hard scale in the 
process and we have to deal with an asymmetric situation. Thus the scales 
so and s' should be chosen as kj alone. At LO accuracy sq is arbitrary and 
we are indeed free to make this choice. Then the arguments of the gluon 
Green's functions can be written as 

S -AL = L^ —=^- (3.12) 

s a a s Pb 

The description in terms of these longitudinal components is particularly 
useful if one is interested in jet production in a hadronic environment. Here 
one can introduce the concept of an unintegrated gluon density in the hadron. 
This represents the probability of resolving a gluon carrying a longitudinal 
momentum fraction x from the incoming hadron, and with a certain trans- 
verse momentum kx- With the help of Eq. (|3.12p a LO unintegrated gluon 
distribution g can be defined from Eq. (|3.4p as 

r-%*p(q) / —x-»Uq,k). (3.13) 



Then we can rewrite Eq. (|3.4p as 
da 



d 2 kjdyj 



d q a / dxi J d q b J dx 2 

x g(x 1 ,<i a )g(x2,q.b)V(qa,xi,q b ,X2;kj,yj), (3.14) 



with the LO jet vertex for the asymmetric situation being 
V(qa,xi,q fe ,x 2 ;kj,yj) = K^ 3 ' (q a , -q 6 ) 

x 5^ (q a + q 6 - kj) 5 L x - ^5e«" j 5 (x 2 - \J^e~ Vj ) • (3.15) 
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Having presented our framework for the LO case, in both 7*7* and hadron- 
hadron collisions, we now proceed to explain in detail what corrections are 
needed to define our cross sections at NLO. Special attention should be put 
into the treatment of those scales with do not enter the LO discussion, but 
are crucial at higher orders. 

3.2 Inclusive jet production at NLO 

A similar approach to that shown in section [3~T1 remains valid when jet pro- 
duction is considered at NLO. The crucial step in this direction is to modify 
the LO jet vertex of Eq. (13. 9p and Eq. (|3.15p to include new configurations 
present at NLO. We show how this is done in the following first subsec- 
tion. In the second subsection we implement this vertex in the symmetric 
7*7* case, and we repeat the steps from Eq. (13, ip to Eq. (|3.1ip . carefully 
describing the choice of energy scale at each of the subchannels. In the 
third subsection hadron-hadron scattering is taken into consideration, and 
we extend the concept of unintegrated gluon density of Eq. (|3.13p to NLO 
accuracy. Most importantly, it is shown that a correct choice of intermediate 
energy scales in this case implies a modification of the impact factors, the 
jet vertex, and the evolution kernel. 

3.2.1 The NLO jet vertex 

For those parts of the NLO kernel responsible for one gluon production we 
proceed in exactly the same way as at LO. The treatment of those terms 
related to two particle production is more complicated since for them it is 
necessary to introduce a jet algorithm. In general terms, if the two emissions 
generated by the kernel are nearby in phase space they will be considered as 
one single jet, otherwise one of them will be identified as the jet whereas the 
other will be absorbed as an untagged inclusive contribution. Hadronization 
effects in the final state are neglected and we simply define a cone of radius 
Ro in the rapidity-azimuthal angle space such that two particles form a 
single jet if R12 = a/ ((f>i - fc) 2 + (yi - IJ2) 2 < Ro- As long as only two 
emissions are involved this is equivalent to the feT _ clustering algorithm. 

To introduce the jet definition in the 2^2 components of the kernel it 
is convenient to start by considering the gluon and quark matrix elements 
together: 



with A2P being the two particle production amplitudes of which only the 
gluonic one also contributes to MRK. This is why a step function is needed 




(3.16) 
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to separate it from MRK. Momentum conservation implies that ki = q a + 

q&-k 2 . 

Expression (|3.16p is not complete as it stands since we should also include 
the MRK contribution as it was previously done in Eq. (12,301) : 

(Kgg+Kqq) (q a ,-q&) = J d D ~ 2 k 2 J dy 2 |#(q a , qt, ki, k 2 )| 2 



J d D 2 k 2 J dy 2 | |A> 9 (q a ,qfe,ki,k 2 )| 2 6>(.SA - s 12 ) 

- /C( B )(q a ,q a - k 1 )/C( B )(q a - ki.-Cfc) ^ (in f| - y 2 ) 

+ |^l 2(? (q a ,q b ,ki,k 2 )| 2 L 



2/2- In ^) 
(3.17) 



We are now ready to introduce the jet definition for the double emissions. 
The NLO versions of Eq. (13. 9p and Eq. (I3.15P then read, respectively, 



[y] 

(a) 



V(q a ,q fe ,7 ? ;k J ,y J ) = [K^ + K^) (q a ,-q 6 ) 

+ y d D ' 2 k 2 dy 2 \B(ci a ,c lb ,kj-k 2 ,k 2 )\ 2 e(Ro-R 1 2) 

+ 2 / d D - 2 k 2 dy 2 |i3(q a ,q b ,kj,k 2 )| 2 0(i? J2 -^o) ^, (3.18) 
J (c) 

V(q a ,xi,q 6 ,x 2 ;kj,yj) = (/C^ +/C^ ) ) (q a , -q 6 ) 



(&) 



fJ 
(a) 



+ / d D ~ 2 k 2 dy 2 \B(c la ,q b ,kj-k 2 ,k 2 )\ z e(Ro-R 12 ) 



(b) 



+ 2 I d 1J - z k 2 dy 2 \B(c la ,c lb ,kj,k 2 )\ z 9(Rj 2 -R ) . (3.19) 



In these two 




L2/J 

(a,6) 




[y] 




(c) 




[x] 




(a,&) " 




[x] 




(c) 



5 (2) (Qa + Q6 - kj) - y^) 
S {2) (q a + q - kj - k 2 ) 5 



v - y {c) ) , 



5 (2) (q a + q - k j) 5 ( x 1 - xi ,b) ) 5 (x 2 -x 



>,&) 



(3.20) 
(3.21) 
(3.22) 
(3.23) 



The various jet configurations demand different y and x configurations. 
These are related to the properties of the produced jet in different ways 



5 (2) (qa + q& - kj - k 2 ) S (xi - x^) 5 ( x 2 - x 
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depending on the origin of the jet: if only one gluon was produced in MRK 
this corresponds to the configuration (a) in the table below, if two particles 
in QMRK form a jet then we have the case (b), and finally case (c) if the 
jet is produced out of one of the partons in QMRK. The factor of 2 in the 
last term of Eq. (13.180 and Eq. (13.190 accounts for the possibility that either 
emitted particle can form the jet. Just by kinematics we get the explicit 
expressions for the different x configurations listed in the following table: 



JET 


y confix 


durations 




x confij 


durations 






a) {- 

b) |<] 




= VJ 

= yj 


x? = 




(a) 
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x 2 ~~ 




-yj 
-yj 


c ) \c 


y (0 = 


ii 4 C) 

x 2 




e yj + lM e V2 

V s 


T ( c ) _ IkjL 


-i/j + 





The variable £ is defined below in Eq. (|3.3ip . Due to the analogue treatment 
of the emission vertex either expressed in terms of rapidities or longitudi- 
nal momentum fractions in the remaining of this section we will imply the 
same analysis for both. In particular, we will not explicitly mention these 
arguments when we come to Eqs. (|3. 411 [3^2]) . 

The introduction of the jet definition divides the phase space into dif- 
ferent sectors. It is now needed to show that the final result is indeed free 
of any infrared divergences. In the following we proceed to independently 
calculate several contributions to the kernel to be able, in this way, to study 
its singularity structure. 

The NLO virtual correction to the one-gluon emission kernel, IC^ V \ was 
originally calculated in Ref . [39l HQl [H] . Its expression reads 



4 V) (q», 



-q&j 



vr^lXl - e) A 2 
3A 2 



(-? + T- 2eC(3) , 



A) i 

+ TT- + 



A 2 



N c e q; 
A 2 

q 2 -q 2 

q b ) 2 + 



-Mi 



hr f 4 

,q& 



A 2 (q 2 + q 2 -4q a q b ) \ /q| 
3(q 2 -q 2 ) 2 J U 2 



A 4 (q 2 +q 2 ) ,j 



6q^q 2 (q^-q 2 ) 2 



(<c + q& -4q a q6) 



(3.24) 



with fa = (11N C - 2n/)/3, C(«) = EfcLi and A = <Ja + 06- ^ £^ 
be expressed in terms of the renormalized coupling constant g^ in the MS 
renormalization scheme by the relation <? 2 = g 2 A r c T(l — e) (47r)~ 2_<E . Note 
that the expression for the virtual contribution given in 95J lacks the log 
squared. 
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Figure 3.3: Feynman diagrams for the process q + q — > q + q + q + q. 



Those pieces related to two-gluon production in QMRK can be rewritten 
in terms of their corresponding matrix elements as 



^QMRK(qa,-qb) = J d D 2 k 2 J dy 2 |^2 9 (qa,q6,ki,k 2 )| 2 6 , (sA - s 12 ) 
gjj^Nl f d D ~ 2 k 2 f 



i 2 ,u- 2e N?_ r d D - 2 

i^b 

and those related to quark-antiquark production are 

|2 



£Q0(qa,-q&) = j d D 2 k 2 J dy 2 |^ 2g (qa,q6,ki,k 2 ) 

gln~ 2e N 2 f d D - 2 k 2 f 
= vr(2vr)^+iq2 q 2 J ^ e(27r) D-4 J ^ Quarto- (3-26) 

For our purposes we need the corresponding amplitudes including those 
parts which vanish in the fully inclusive case. Therefore we have calculated 
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these amplitudes at which the ingoing legs are off-shell. To obtain them in 
a gauge invariant way it is necessary to consider an appropriate embedding 
as a quark-quark scattering where the gluon or quark-antiquark pair is 
produced additionally. The according Feynman diagrams are depicted in 
Figs. 13.21 and 13.31 and using the Mandelstam invariants s, t, and u, and we 
extracted the following results 
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These expressions are in agreement with the corresponding ones obtained in 



44 



Chapter 3. Jet production 



Ref. [95]. The following notation has been used: 
|ki| 



(3.29) 



|ki| + |k 2 | e ^y ' 

A =(1 — x)k\ — ric 2 , (3.30) 
A 2 

S =s + A 2 = — + A 2 , (3.31) 

x(l — x) 

1 — 2r 

E= 2(2x - l)q 2 + 4Aq a + " 



x(l — x) 



n 2 

2 i o * A \ 



= -2x(l-x)((2x-l)A 2 + 2AA) — ^ la- ( 3 - 32 ) 

v 7 x(l - x)A 2 + A 2 

We now study those terms which contribute to generate soft and collinear 
divergences after integration over the two-particle phase space. They should 
be able to cancel the e poles of the virtual contributions in Eq. (|3.24|) . i.e. 



p=4„-2e A 

r {v) , \ 9^ 4_ 

^singular IVtaj Q6j — _ A ^2 




(3.33) 



Here we identify those pieces responsible for the generation of these poles. 

One of the divergent regions is defined by the two emissions with mo- 
menta k\ = a±pA + PiPb + k\± and fc 2 = a 2 pA + P2PB + k 2 ± becom- 
ing collinear. This means that, for a real parameter A, k\ ~ A&2, %■£■ 
ki± ~ A/c2±, Oi\ — Aa2 and thus a2&i_i_ — aik 2 ± — 0. Since x = a 4^, 
this is equivalent to the condition A ~ 0. In the collinear region s = x ^_ x ^ 
tends to zero and the dominant contributions which are purely collinear are 
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(3.35) 



The quark-antiquark production does not generate divergences when ki or 
k2 become soft, therefore we have that the only purely soft divergence is 



a singula] 
gluons 



+ W(3)+%) - ^ (3-36) 



soft " a 10 V Atii kf k 2 4kf k 



2 



where we have used the property that, in the soft limit, the tu product tends 
to q 2 q 2 . We will see that these terms will be responsible for simple poles 
in e. The double poles will be generated by the regions with simultaneous 
soft and collinear divergences. They are only present in the gluon-gluon 
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production case and can be written as 
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Focusing on the divergent structure it turns out that in the soft and collinear 
region the first line of Eq. (|3.37p . ^4(5), has exactly the same limit as the 
second line, . This is very convenient since we can then simply write 
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The MRK contribution of Eq. ([XT7]) has the form A M rk = - 
when added to all the other singular terms we get the expression 



4^4(4) and 



dy 2 |£ s (q a ,qfe,k 2 ,ki 
d D - 2 



^(2vr) 



"^2 f r f a singular a , \ . singular 1 

^4 J {^gluons ( fi A - S12J + ^ qU a rks j, 

(3.39) 
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(3.40) 



We have labeled the different terms to study how each of them produces the 
e poles. We will do this in section 13.31 

With the singularity structure well identified we now return to Eqs. (|3.18l 
I3.19P and show that they are free of any divergences. Only if the di- 
vergent terms belong to the same configuration, this cancellation can be 
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shown analytically. With this in mind we add the singular parts of the 
two particle production of Eq. (|3.39|) in the configuration (a) multiplied by 
= 1- 9(R - R l2 ) - 9(R 12 - Rq): 



V 



fc B )+/CM N )(q a ,-q & ) 



+ J d D 2 k 2 dy 2 |£ s (q a ,q b ,k/ - k 2 ,k 2 
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9(Rq — Ri 2 ) 



(a) 



. (3.41) 



The cancellation of divergences within the first two lines is now the same 
as in the calculation of the full NLO kernel. In section [3731 we will show how 
the first two lines of Eq. (I3.4ip are free of any singularities in the form of 
e poles. In doing so we will go into the details of the role of s\. The third 
and fourth lines are also explicitly free of divergences since these have been 
subtracted out. The sixth line has a ki <-> k 2 symmetry which allows us to 
write 
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(c) 



|^(q a ,q fc ,kj - k 2 ,k 2 )| 2 0(i? 12 - flo)0(|ki| 



. (3.42) 



We can now see that the remaining possible divergent regions of the last line 
are regulated by the cone radius Ro. 

It is worth noting that, apart from an overall a 2 (n 2 ) factor, the NLO 
terms in the last four lines in Eq. (|3.42[) do not carry any renormalization 
scale dependence since they are finite when e is set to zero. The situation 
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is different for the first two lines since V contains a logarithm of fi 2 in the 
form 

V = V™ ( 1 - In 41 + AV. (3.43) 

where AV contains the third to sixth lines and the /^-independent part of 
the first two lines of (|3.42p . It is then natural to absorb this term in a 
redefinition of the running of the coupling and replace a s (fj?) by a s (k 2 ). 
For a explicit derivation of this term we refer the reader to section 13.31 

Therefore we now have a finite expression for the jet vertex suitable for 
numerical integration. This numerical analysis will be performed elsewhere 
since here we are mainly concerned with the formal introduction of the jet 
definition and the correct separation of the different contributions to the 
kernel. 

What remains to be proven is the cancellation of divergences between 
Eq. (|3,33[) and Eq. (|3.39|) . This will be performed in section [3~B1 Before doing 
so, in the next two subsections, we indicate how to introduce our vertex in 
the definition of the differential cross section. Special care must be taken in 
the treatment of the energy scale in the reggeized gluon propagators since in 
the symmetric case it is directly related to the rapidity difference between 
subsequent emissions, as we will show in the next subsection, but in the 
asymmetric case of hadron-hadron collisions it depends on the longitudinal 
momentum fractions of the t-channel Reggeons. 



3.2.2 Production of jets in 7*7* scattering 

We now have all the ingredients required to describe the inclusive single 
jet production in a symmetric process at NLO. To be definite, we consider 
7*7* scattering with the virtualities of the two photons being large and of 
the same order. All we need is to take Eq. (|3.4p for the differential cross 
section as a function of the transverse momentum and rapidity of the jet. 
The vertex V to be used is that of Eq. (|3.42p in the representation based on 
rapidity variables of Eq. (|3.18p . The rapidities of the emitted particles are 
the natural variables to characterize the partonic evolution and s-channel 
production since we assume that all transverse momenta are of the same 
order. 

Let us note that the rapidity difference between two emissions can be 
written as 

m - y i+1 = In (3.44) 

which supports the choice SR-i^+i = ^k 2 \s. 2 +1 in Eq. (|2.15p . This is also 
technically more convenient since it simplifies the final expression for the 
cross section in Eq. (|2.25p . 
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In Fig. l3.4l we illustrate the different scales participating in the scattering 
and the variables of the evolution. We give the conditions for MRK: all 
transverse momenta are of similar size and much larger than the confining 
scale, the rapidities are strongly ordered in the evolution from one external 
particle to the other. At each stage of the evolution the propagation of the 
Reggeized gluons, which generates rapidity gaps, takes place between two 
real emissions. There are many configurations contributing to the differential 
cross section, each of them with a different weight. Eq. (|3.4p represents the 
sum of these production processes. 
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Figure 3.4: Momenta for 2 — ► 2 + (n — 1) + jet amplitude in the symmetric 
configuration with MRK. The produced jet has rapidity y,j = yj 



and transverse momentum kj = k,-. 



3.2.3 The unintegrated gluon density and jet production in 
hadron— hadron collisions 

In this subsection we now turn to the case of hadron collisions where MRK 
has to be necessarily modified to include some evolution in the transverse 
momenta, since the momentum of the jet will be much larger than the typical 
transverse scale associated to the hadron. 

In the LO case we have already explained that, in order to move from 
the symmetric case to the asymmetric one, it is needed to change the en- 
ergy scale from the choice in Eq. (13.70 to the one in Eq. (13.120 . This is 
equivalent to changing the description of the evolution in terms of rapid- 
ity differences between emissions to longitudinal momentum fractions of the 
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fc-ordered MRK: 



p\ < kf < . . . < k|_! < k] (Jet) k| > kj +1 • • ■ > k^ > p| ~ A 

Va = V0>Vi> ■■> Vj-i > Vj > Vj+i ■>Vn> Vn+l = 1 
To impact factor 

* To kernel 

V - \ 



2 

QCD 




q l ' up ^ k l ' down I / 2 \ j 



(V1-V2) fsv (say 



(I) 1 



7 up \ * down 



To emission vertex 



Vq„ / down V k n-1 



k 2 



+ 1 / down 



Figure 3.5: Momenta for 2 — > 2 + (n — 1) + jet amplitude in the asymmetric 
configuration with /^-ordered MRK. 



Reggeized gluons in the i-channel. Whereas in LO this change of scales has 
no consequences, in NLO accuracy it leads to modifications, not only of the 
jet emission vertex, but also of the evolution kernels above and below the jet 
vertex. These new definitions will allow the cross section still to be written 
in a factorizable way and the evolution of the gluon Green's function still to 
be described by an integral equation. 

To understand this in detail we start by writing the solution to the NLO 
BFKL equation iteratively, i.e. 



f 1 °° P _1 f 1 



(3.45) 



where qi = k a and qj = q a . As both sides of the evolution are similar, 
we now focus on one side of the evolution towards the hard scale and use 
Fig. 13.51 as a graphical reference. Starting with the symmetric case, the 
differential cross section for jet production contains the following evolution 
between particle A and the jet: 



da 



d 2 kjdyj 



« 2 - / 



x /^ / " (k »' q ° ) 



SAJ 



V(q a ,q 6 ;kj,yj)... (3.46) 
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In the asymmetric situation where k 2 3> k 2 the scale wk 2 kj should 



be replaced by k 2 . In order to do so we rewrite the term related to the 
choice of energy scale. To be consistent with Fig. 13.51 we take kj = kj, 
ko = — k a = — qi and qj = q a . To start, it is convenient to introduce a 
chain of scale changes in every kernel: 



saj 



n 

i=l v 



(3.47) 



which can alternatively be written in terms of the t-channel momenta as 

ri- 1 



SAJ 
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,i=i 



n 2 \ 2 



q« 



(3.48) 



For completeness we note that we are indeed changing the variable of evo- 
lution from a difference in rapidity: 



saj 



(3.49) 



to the inverse of the longitudinal momentum fraction, i.e. 



saj __ _J_ 
k 2 a/ 



(3.50) 



This shift in scales translates into the following expression for the cross 
section: 
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(3.51) 



As we mentioned above these changes can be absorbed at NLO in the kernels 
and impact factors, we just need to perturbatively expand the integrand. 
The impact factors get one single contribution, as can be seen in Fig. 13.51 
and they explicitly change as 

£(k ) = <K(k a ) - ~k* J d 2 q^^/C^(q,k a )ln^. (3.52) 
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The kernels in the evolution receive a double contribution from the differ- 
ent energy scale choices of both the incoming and outgoing Reggeons (see 
Fig. I3.5p . This amounts to the following correction: 



/C( qi ,q 2 )=/C( qi ,q 2 )-i / d 2 q/C( B )( qi ,q)/C^(q,q 2 )ln4- (3-53) 

2 J q 2 

There is a different type of term in the case of the emission vertex where the 
jet is defined. This correction has also two contributions originating from 
the two different evolution chains from the hadrons A and B. Its expression 
is 

2 

2 



V(q a ,q fe )=V(q a ,q b )-i / d 2 q)C^ (q a , q}V^ (q, q&) In - — — 

- \ [ d 2 qV^(q a ,q)/C( B )(q,q b )ln q2 (3.54) 
2 J (q a - q) 2 

These are all the modifications we need to write down our differential 
cross section for the asymmetric case. The final expression is 

d^jdyj-J dCla J dK Jfcrkg 

As in the LO case, we can use Eq. (|3.50p to define the NLO unintegrated 
gluon density as 

g(x,k) = J d\^0L J ^(k.q)^. (3.56) 

The gluon Green's function f u is the solution to a new BFKL equation with 
the modified kernel of Eq. (|3.53p which includes the energy shift at NLO, 
i.e. 

w/ w (k a ,q a ) 

= <y (2) (ka - q«) + j d 2 q/C(k a ,q)/ w (q,q a ). (3.57) 



In this way the unintegrated gluon distribution follows the evolution equa- 
tion 

dg(x,c[ a ) 
dlnl/x 

Finally, taking into account the evolution from the other hadron, the differ- 
ential cross section reads 

da 



yd 2 q/C(q a ,q) 5 (x,q). (3.58) 



d 2 kjdyj 



J d 2 q a J d 2 q b g(x a ,ci a ) g(x b ,ci b ) V(q a , q fe ; kj, yj), (3.59) 
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with the emission vertex taken from Eq. (|3.54p . 

We would like to indicate that with the prescription derived in this sub- 
section we managed to express the new kernels, emission vertex and impact 
factors as functions of their incoming momenta only. It is also worth men- 
tioning that the proton impact factor contains non-perturbative physics 
which can only be modeled by, e.g. 



where pi are positive free parameters, with Qq representing a momentum 
scale of the order of the confinement scale. Though Regge theory suggests 
a P2 ~ 0.08 [96] the common collinear parton distributions available work 
with a rather valence like value of p2 ~ —0.09 |97j . But in any case the 
initial x dependence in this expression would be of non-perturbative origin. 

Let us also point out that the prescription to modify the kernel as in 
Eq. (|3.53p was originally suggested in Ref. [18] in the context of deep inelastic 
scattering. This new kernel can be considered as the first term in an all 
orders perturbative expansion due to the change of scale. When all terms are 
included, the kernel acquires improved convergence properties and matches 
collinear evolution as shown in section [2.3.2L 

3.3 Cancellation of divergences and a closer look 
at the separation between MRK and QMRK 

During the calculation of a NLO BFKL cross section, both at a fully inclusive 
level and at a more exclusive one, there is a need to separate the contribu- 
tions from MRK and QMRK. In order to do so we followed Ref. |37| and 
introduced the parameter s\ in Eq. (I2.20P and Eq. (|2.2ip . In principle, at 
NLO accuracy, our final results should not depend on this extra scale. In 
fact, as we have remarked earlier in our discussion of the total cross section 
(after Eq. (|2.30p ). we could have taken the limit s\ — > oo: the logarithms 
of sa cancel, and the corrections to the finite pieces die away as ©(s^ 1 ). In 
the context of the inclusive cross section, however, we prefer to treat sa as 
a physical parameter: it separates MRK from QMRK and, hence, cannot 
be arbitrarily large. We will therefore retain the dependence upon sa: in 
the remainder of this section we demonstrate that, in our inclusive cross 
section, all logarithmic terms cancel (analogous to Eq. (|2.30p ). and we leave 
the study of the corrections of the order O (s^ 1 ) for a numerical analysis. 
It will also be interesting to see how this dependence on sa could be related 
to the rapidity veto introduced in Ref. [551 1591 l9"51 19"9] . 

Let us consider the sa dependent terms in Eq. (I3.39P which are only 




(3.60) 
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present in the gluon piece: 
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IV 



0(s A -s 12 ) = J^cSi, (3.61) 



where we have used the numbering (J, II, IV) corresponding to, respec- 
tively, (Gluon| co ii 1 ,Gluon| co ii 2 ,Gluon| so f t ,Gluon| soft&co ii) in Eq. (|3.40l) . 

To calculate each of the Si terms we start by transforming the rapidity 
integral into an integral over x in the form J dAy = J x ^ x ^ ■ We consider 
s\ much larger than any of the typical transverse momenta. In the limit of 
large sa, the theta function 6(sa — s) amounts to the limits j 1 + O (s^ 2 ) 

and 1 — ^ + O (s^ 2 ) for the x integral. 
We firstly consider Sni which is 
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We are only interested in the logarithmic dependence on s\ and hence we 
do not need to calculate O (s^ 1 ) or s\ independent factors. 
The next s\ contribution we calculate is Sjy which reads 
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(3.63) 



The part with logarithmic s\ dependence can be calculated analytically: 
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It is then clear that this logarithmic sa contribution cancels against that of 
S IU in Eq. 453511 ■ 

Let us proceed now to show that the contribution of Si is directly of 
O (s^ 1 ) and does not contribute with any logarithm of s\. In the relevant 
integral we introduce the change of variable k 2 — > A = (1 — x)A — k 2 and 
obtain 
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(3.65) 



We do not write here the lengthier, but similar expression which corresponds 
to Si 1 and also only contributes to O (s^ 1 ). 

With this we have shown that the sum of different terms in Eq. (|3.61|) 
is free of logarithmic dependences on sa proving, in this way, that the re- 
maining O (s^ 1 ) corrections vanish at large values of s\. In particular, it is 
possible to take the s\ — » 00 limit in order to completely eliminate the de- 
pendence on this scale. This is convenient in the fully inclusive case where it 
is very useful to write a Mellin transform in the kx dependence of the NLO 
BFKL kernel. 

If we perform this s\ —* 00 limit, then Sin and Siy can be combined 
and their sum is 
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(3.66) 
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Regarding Si, one obtains from the integration 
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The contribution from Sn is more complicated and the relevant integral can 
be obtained in the following way: 
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We now need to integrate it over x to obtain: 

1 vr vr(l + e)r(2 + e) 



Sn 



A 2 

M 2 



(3.68) 



(3.69) 



A 2 (4tt) £ sin(7re)r(4 + 2e) 

This result gives the same poles in e as the result given in [95], but differs 
for the finite contribution. To obtain all the e poles we now also include the 
quark contributions present in Eq. (|3,39p . We denote them as 

• singular VI 



d D ~ 2 k 2 

^(2tt) d -^ 



dy2 



A 



quarks 

q 2 q& 



i=V 



(3.70) 



where the correspondence with Eq. (|3.40p is (V, VI) — > (Quark^d^ , Quark| co n 2 ). 
Adding everything up, the sum of all the terms reads 



VI 



1 TiT(l-e) 



A 2 



(47T) 
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M 2 
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2N c e 
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The final expression for Eq. f|3.39[) is thus given by 



d D 2 k 2 / dy 2 |S s (q a ,q6,kj 



ko, ko 



k 2 
K ./ 



A ^1 10n f 



(3.72) 



When we combine this result with the singular terms of Eq. (13.240 then 
we explicitly prove the cancellation of any singularity in our subtraction 
procedure to introduce the jet definition. The finite remainder reads 



a?V) 1 



7T 



k 2 



+ 



1 

12 



2tt z 



■ Po 



(3.73) 



We have already discussed the logarithmic term due to the running of the 
coupling in Eq. (|3.43p . The non-logarithmic part is similar to that present in 
other calculations involving soft gluon resummations [100, 1101] where terms 
of the form 

a s (l + Sa s ) (3.74) 

appear and offer the possibility to change from the MS renormalization 
scheme to the so-called gluon-bremsstrahlung (GB) scheme by shifting the 
position of the Landau pole, i.e. 



Agb = AMs ex P S 



.2N r 



The factor S differs from ours in the 7r 2 term: 



1 



12 



0o 



S = — 4 - n z + 5^ . 



(3.75) 



(3.76) 



The origin of this discrepancy lies in the fact that we used the simplest form 
of subtraction procedure. In appendix lAl we suggest a different subtraction 
term which is more complicated in the sense that it subtracts a larger portion 
of the matrix element in addition to the infrared divergent pieces. When 
this is done and we put together the divergent pieces of Eq. (|3.24|) and the 
second line of Eq. (|A.12p , then we recover the same S term. 



Chapter 4 

Angular decorrelation 



After the construction of a jet vertex in the central region at NLO, we 
dedicate this chapter to another aspect of jet physics, namely its angular 
correlation. 

In the first section, we start with a short review of the LO description of 
Mueller-Navelet jets and introduce some notation, while in the second sub- 
section the extension to NLO is presented. We then study the phenomenol- 
ogy and compare our results with experiment. In this context, we discuss 
the need for a partial resummation of even higher order contributions and 
show its consequences. Thereafter, we turn to the case of electron-proton 
scattering. The angular decorrelation between the electron and the most 
forward jet is governed by the same physics as in the Mueller-Navelet case. 
To make use of our experiences from proton-proton collisions for a study 
of angular correlation in deep inelastic scattering (DIS) is hence a natural 
extension which will be presented in the second section of this chapter. 

4.1 Angular correlation of Mueller-Navelet jets 

A discussion of the Mueller-Navelet jet has been already given in the context 
of the power like rise of the cross section in section [2^1 Our focus, now, lies 
on the angular correlation of the jets as a more exclusive observable. Before 
we go into the details, it is worthwhile to point out that the Mueller-Navelet 
jets lie at the interface of collinear factorization and BFKL dynamics. The 
partons emitted from the hadrons carry large longitudinal momentum frac- 
tions. After scattering off each other, they produce the Mueller-Navelet jets. 
Because of the large transverse momentum of the jets, the partons are hard 
and obey the collinear factorization. In particular, its scale dependence is 
given by the DGLAP evolution equations. Between the jets, on the other 
hand, we require a large rapidity gap which is described by BFKL dynamics. 
The hadronic cross section, hence, factorizes into two conventional collinear 
parton distribution functions convoluted with the partonic cross section, to 
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be described by the BFKL equation. With respect to the partonic cross 
section, the incoming partons, consequently, are considered to be on-shell 
and collinear to the incident hadrons. 

For the angular correlation, theoretical predictions from LO BFKL exist 
|214l22j. and even improvements due the running of the coupling and proper 
treatment of the kinematics have been implemented [231 [2"4"j . A first step 
to a NLO description has been made in Ref. |32} [33] on which our work - 
starting with the second subsection - builds up. But for pedagogical reasons 
we start with the LO description. 



4.1.1 Mueller-Navelet jets at LO 

Our starting point is the total cross section of two particles as given in 
Eq. (|3.ip . For the case of Mueller-Navelet jets these particles are partons 
out of the protons. We initially consider the partonic gluon-gluon scattering: 

f cPqi f d 2 q 2 [did , w , N wAv 

a{s) = J 2Wj 2^q|7 ^-^on,i(qi)^iuon, 2 (q 2 )/.(qi,q 2 )e ». 

(4.1) 

In this symmetric situation we have chosen so such that it leads to a de- 
scription in terms of the rapidity Ay = In | qi | S | q2 | between the two scattered 
gluons. To be observed by the detector as a jet, a gluon has to have a min- 
imal transverse momentum. We incorporate this resolution scale p 2 by a 
Heaviside function defining the jet vertex: 

«W(qi) - ^et,l(qi) = %f^^(q?-p?). (4.2) 

Simultaneously we rearranged the constants in a convenient way and, again, 
make use of the convention a s = a s N c /ir. At LO we are completely free in 
our choice of sq. The introduction of the resolution scale for the jets suggests 
a more convenient rapidity variable Y = In — ^— . We will comment on this 
change in more detail at the next section dealing with the NLO description. 

According to the normalization of the BFKL equation used in this work, 
the gluon impact factor at LO in our case reads 



TTOt 

*gluon,i(qi) = (4.3) 



Due to high-energy factorization the differences between quark-quark, quark- 
gluon, and gluon-gluon scattering concern only a constant factor in front of 
the impact factors, namely 

Cp 

$quark,l(qi) = ^*tfuon,l(qi), ( 44 ) 

where Ca = N c and Cf = (N 2 — 1)/(2N C ) are the Casimir invariants of 
the adjoint and fundamental representation, respectively. Hence, we can 
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incorporate the combined quark and gluon contributions to the hadronic 
cross section by the effective parton density [102 

/eff0c,/4) = G(x,n 2 F ) + \ V [Qf{x,n\) + Qf(x, fJ 2 F )] , (4.5) 

9 / 

where the sum runs over all quark flavors, and \ip denotes the factorization 
scale. For the proton-proton scattering the Mueller-Navelet cross section 
then reads 



dx\ dx2 



f e s(xi, n F )f eS (x 2 , n F )a(xix 2 s). (4.6) 



For the following steps we focus on the partonic cross section, and the 
convolution with the parton densities will be performed at the end. The 
kinematics are depicted in Fig. 14.11 and we can write the partonic cross 
section (|4.ip as 



Hs) J d 2 ki J d 2 k 2 J ^(k 1 |l» jct , 1 / aJ |. jet , 2 |k 2 )e- y (4.7) 

f; J d 2 k! J d 2 k 2 J^lJ d\, J d\ 2 JdvJ dv> 



2 

n,n'=— oo 



(ki 1 4>jet,i | qi) (qi | v, n) (n, v \ f w \ v' , n') (n' , v' \ q 2 ) (q 2 1 *j et ,2 1 k 2 ) e wY . 

(4.8) 



Here we made use of the operator representation introduced in section 12.31 
The integrals over transverse momenta have to be taken over the whole 
two dimensional space while the v integrations go from — oo to oo. The 
contour of the U) integration is to the right of all w-plane singularities of the 
integrand. 

Due to the cylindric symmetry around the beam axis, the cross section 
can only depend on the difference between the azimuthal angles of k\ and 
k 2 , but not on their absolute values. Nevertheless, it is convenient for our 
derivation to fix an arbitrary axis in the azimuthal plane and to denote the 
azimuthal angles of k\ and k 2 with a\ and a 2 , respectively. In the end, we 
are interested in the angle <f) = a% — a 2 — it between both jets. Therefore, 
we transform the integral measure in the following way and later integrate 
one of the angles (a 2 ) out (see Eq. (14.13H ) to give a factor of 2ir. 

d 2 kl d 2 k 2 d\, d 2 q 2 = da2 2 dk * dW d\ 2 

= de P da 2 d ^ d2 * 1 ' k2 ' 2q2 . (4.9) 



The integrations belonging to the second jet vertex can be easily per- 
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Figure 4.1: Kinematics of the partonic cross section in the Mueller-Navelet 

case 



formed. The result leads to the following definition of the coefficient C2{y'\- 



dk 2 / drqz (n ,v |q2> (q 2 |*j etj2 |k 2 ) 
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2 1 dkl j d2 ^,^2 



/12 
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k 2 
K 2 
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s/2\ + 



IV 



! -in a 2 



2tt 



=: <*(!/)- 



-471 Q2 



2tt 



(4.10) 



When performing the same calculation for the first jet vertex, we have to 
note that, due to conventions, qi is an incoming momentum with respect to 
the jet vertex while qi is an outgoing one. This means that qi has to equal 
— ki, and hence we get for the first jet vertex 



- I dk( I cPqi (ki|$ jct ,i|qi)(qi|i/,n) 



1 1 



\fl\-iv 



\Pi) — — = : c n^)- 



2vr 



2vr 



Recalling Eq. (|2.4ip we have 



dio 



{n,v\f w \v ,n)e' 



2iri u - a s xo [\n\,\ + iv) 

a sX o(\n\ ,\+iv)Y S(v-i/)6 nn .. 



(4.11) 

5{v - v')8 nn ,e» Y 
(4.12) 
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Inserting these results in Eq. (|4.7p leads to the following compact ex- 
pression 



2-2 00 
7Ta« 
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n,n — — oo 



da 2 / du I du' 



d(u 



s in(ai— 7r) 
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2tt 



9 _ 9 OO 

7ra; 
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Ei V" p a«Xo(|n|,|+^)^ 



from which we can read off the angular distribution 
<t(Y) ir 2 a 2 s 



4 V^Mn^oo 



with 



C^ L (E) 



2vr 1 + ^ ^2 



dt/ 1 ( ri\ V p a sX0 {\n\,±+i»)Y 



(4.13) 



(4.14) 



(4.15) 



This decomposition is very well suited to study the angular correlation. 
Besides the differential cross section itself, the mean values < cos(n</>) > can 
be easily obtained from Eq. (|4.15|) due the orthogonality of the functions 
exp(in^): 

<cosM)>=^|. (4.16) 



4.1.2 Mueller-Navelet Jets at NLO 

In general, a complete NLO calculation modifies the jet vertices <3?j e t and the 
Green's function f w as well. The decorrelation between the jets is assumed 
to be mainly caused by the Green's function even if the NLO jet vertex 
contains one additional emission. In particular the rapidity dependence is 
exclusively driven by the Green's function. Therefore, we focus on the effect 
of the NLO Green's function while keeping the jet vertices at LO. 

Nevertheless, we have to state that the NLO Mueller-Navelet jet vertex 
in principle has been calculated [251 [26] . Since it is the link between the 
collinear factorization from the hadron side and the BFKL dynamics be- 
tween the jets, the NLO calculation reveals divergent contributions which 
have to be combined with the NLO collinear parton density and others which 
overlap with the BFKL evolution. It has been shown in Refs. [25} I26j that 
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after all the factorization remains valid with a finite result for the jet ver- 
tex. However, the final results are not yet in a form appropriate for a direct 
imp lenient at ion . 

Already in the previous chapter we stressed that the choice of sq at 
NLO is no longer completely free. In the LO discussion we made use of 
the possibility to perform the qi/ 2 and ki/2 integrations without touching 
the Green's function. This was possible since we switched from the physical 
rapidity Ay to an effective rapidity variable Y. A consistent NLO calculation 
demands to compensate for this energy scale change 



s 



)""(■ 



X\X2S 



^Ay 



by a change of the impact factors 



s \ u _ ( XlX 2 S \ 



■jY 



(4.17) 



2 J Pi 

(4.18) 

The superscript (Ay) marks the correct expression for the energy scale sq 
and (Y) that for the energy scale s' . The Born expressions are labeled with 
a B. Therefore, this scale change is a NLO effect to the impact factors which 
we intend to postpone to a future study. As a matter of fact, this choice of 
energy scale separates the integration over the transverse momenta from that 
over the longitudinal momentum fractions and from the rapidity itself. Due 
to this separation the integration over the longitudinal momentum fractions 
factorizes and just leads to an overall factor. 

To calculate (n, v\ fj\v' , n'), we first have to determine the matrix ele- 
ments (n,v\fc\v',n') with the help of Eq. f|2.42[) : 



{n,v\K,\u ,n ) = / d qi / d c& (n, Hqi)(qi|£|q2)(q2|n v 



\n ( \n'\, - + iv' I + l 



o I In |, - +iv 



(4.19a) 



We can perform the qi integration if we take into account the explicit form 
of the running of the coupling a s (qf) = a s (y?) — a 2 s (y?)-^- In The 
logarithm of q| can be expressed as a derivative with respect to either v or 
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v 1 . We choose a version, symmetrized in v and v': 
(n, u\K\u , n) 



2vr 2 



Xo I \n |, - + iu 



Xo n , r + w 



4N, 



+ — — 5 \\n h^+ w 



. d . d . , 



Xo [ \n |, - + %v I + a s xi I \n |, - + iv 



d d 

+ *a 21n/i z 

ov' ov 



+ ^ 



a s /3o Xo § + ii/) 



8iV, 



(4.19b) 



For the last transformation we used the explicit form of the function 
S(n, 7) given in Eq. (|2.43p and, from now on, a s = a<j(/i 2 ). The derivative 
operators inside the curly brackets now act on the Dirac delta distribution 
which - at the end - has to be resolved by integration by parts. The v 
derivative, with flipped sign, will then act on v dependent terms outside 
this kernel matrix element. The same holds for the v' derivative, but it will 
also act on the leading order kernel attached to this derivative operator. 
That term cancels the derivative of %o from the NLO kernel. This is very 
convenient since it is the only term on the r.h.s. of Eq. (|2.42p which is not 
a real valued function for real v or, in other terms, which is not symmetric 
under the exchange 7 «-> 1 — 7- Since in the full expression of the cross 
section the v dependent part is to the left of the Green's function and the v' 
dependent part to the right, we already evaluate the Dirac and Kronecker 
deltas by explicitly indicating that one derivative is acting to the left and 
the other one to the right. Therefore, we can write the diagonal matrix 
elements of Eq. (|4.19bp without the Dirac and Kronecker deltas as 

a s p 



O.s- 



Xo 



8N r . 



{-idvXo + iXodu - 21n/i 2 | + a s xi 



(4.20) 



where we introduced the short hand notations X0/1 = Xo/i(l n l> + iv). 

As shown in chapter [2j the operator of the gluon Green's function is 
given as the inverse of another operator / w = (u — (see Eq. (I2.36|) ). To 
calculate the matrix element (n, v\f^\v' , n'), we there simply used Eq. (|2.40D 
to write 

L = {u-KL)- 1 => (n,v\fl\v',n') = * (4.21) 

u — {n, v\rL\v , n) 
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One might pose the question whether an operator in the denominator of a 
fraction is well defined. From the mathematical point of view, this is not 
trivial but well defined since that operator inversion actually is given by the 
Neumann series 



, h 



A, = (w-^- 1 = ^E(-] • (4.22) 

Since the kernel acts on its eigenfunctions, the eigenvalues stand in for the 
kernel, and the series becomes a simple geometric series. That series can 
again be written as one simple fraction. Due to the running of the coupling, 
the NLO kernel is not really diagonal with respect to the LO eigenfunctions 
such that the "eigenvalues" are still operators. Therefore, we have to study 
the action of the ^-derivatives in the series representation of (n, v\ fj\v', n'). 
At this stage we can write the partonic cross section as 

n=-oo J J J k=0 V 7 

(4.23) 

where X is defined in Eq. (|4.2U|) . The derivative operators within X now 
act on other X's and on the jet vertex coefficient functions c\{v) and c-iiy} ■ 
For every operator d v we find an operator d v at the mirrored position which 
acts on the exactly mirrored bunch of X's. Since the operators, acting to the 
right, come with a different sign compared to those, acting to the left, these 
contributions cancel. The only exception emerges when the derivatives act 
on the coefficient functions Cxjiiy)- Therefore, the operator d u effectively 

can be replaced by d 2 {v) I 'c-iiy) if the according operator d u at the mirrored 
position is replaced by d^{y)jc\{y) at the same time. In this way we can 
evaluate all derivatives such that every X is replaced by a simple function: 



d lnCl H 



du c 2 (f) 



21n// 2 |>+a s xi- (4.24) 



For the case of our jet vertices defined in Eqs. (14.101 14. lip the logarithmic 
derivative in Eq. (|4.24[) reads 

-4 ln ^TT = ln ^) + T^- (4-25) 

Now the operator series in Eq. (|4.23p reduces to a simple geometric series, 
and it is straight forward to perform the Mellin transformation to calculate 
the NLO coefficient 



-oo 2vr \p 2 2 ) \ + u 2 

(4-26) 
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In the derivation of Eq. (|4.19p we expressed the coupling constant at the 
scale in terms of one at a general renormalization scale fi. Since this was 
an intermediate step, it is natural to redo it for the final result. Thus we 
can reexpress the running coupling term in Eq. (|4.26|) in the following way 

-2 A) , PlP2 - f \ I A Q7\ 

as ~ as W~ ~ 2~ — ' a s{PlP2), (4.27) 
and obtain as coefficient 



e L on = rp(4Y - — rr-r — ( 4 - 28 ) 

In section [2.3.2l we have shown that it is necessary to resum terms of even 
higher order than NLO to obtain a reliable kernel. The resummed kernel, 
including the impact factor contribution — , then replaces the kernel 

in the exponent of Eq. (|4.28|> and will be referred to it as 

r°° rhj fr, 2 \ iV ^ rcsum (n,i/)Y 

CT nm (Y)= r- 2 ■ (4.29) 

The technical details, how to include this additional contribution to the 
kernel in the resummation procedure, can be found in appendix IB.31 

It is not possible to perform the ^-integration in Eq. (|4.29|) analytically. 
In principle there are two ways two evaluate it: either by a numerical in- 
tegration or by an analytical approximation. In our studies we perform 
the integration numerically although this is quite time consuming as the 
resummed kernel is given by an implicit equation. 

The saddle point method is the common way to approximate this kind of 
integrals analytically. It relies on the assumptions that the exponent u(n, v) 
is well approximated by a Taylor expansion around its maximum at v = 

uj(n, v) « u{n, 0) + - \ U=Q , (4.30) 

and the multiplicative function in front of the exponential is 'slowly varying' 
around this maximum. If we write a general kernel as u(n, v) = a s x(n, 1/2+ 
iv) and make use of these approximations, the saddle point approximation 
to Eq. (P~29|) reads 



a a (piP2)Y[ xo+a s (PiP2) Xi-glrTT 1 



C n (Y) « ^ /27T ^ X ^ )Y 
TT\j a s x" {n,^)Y 

where x"{ n il) = d 2 x(n,^i)/d^ 2 . 

For the LO kernel this approximation works quite well, but for the pure 
NLO kernel it breaks down since x"( n A/2) < in that case. Even after 
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a s x{m = 0,7 = \) 




Figure 4.2: The kernel and its second derivative in dependence of the cou- 
pling a s . The solid lines show the results for the LO (straight 
line) and NLO (bended line) kernel without resummation. The 
results for the different resummation schemes lie between these 
both. In both diagrams they are from top to bottom: scheme 4, 
scheme 3, scheme 1, scheme 2. 
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the resummation procedure the saddle point method is not sufficient for 
the resummed kernel. To illustrate this, we show in Fig. 14.21 the resummed 
kernel and its second derivative with respect to 7 = 1/2 + iv. Since the 
resummed kernel is implicitly 00 dependent, we have 



X \n 



uJ 



d 2 x(n^,u) 



di 2 



1 



a, 



dx {n,J,uj) 
doj 



dj 2 



(4.32) 

Fig. 14.21 illustrates that the kernel at the saddle point for the resummed ver- 
sion is positive also for larger values of a s while the NLO kernel without re- 
summation turns and becomes negative. The different resummation schemes 
give compatible results for the kernel. In contrast, the second derivative 
shows a stronger dependence on the resummation scheme. In scheme 1 and 
2 the maximum even turns into a minimum for larger a s as it is the case 
for the NLO kernel without resummation at a s ~ .05. Already at a s = .1 
the second derivatives deviate from each other significantly at the saddle 
point, but the complete kernel is very similar for the different resummation 
schemes, as it can be seen in Fig. 14.31 Therefore, we conclude that it is 
mandatory to perform the ^-integration numerically to obtain reliable re- 
sults. As Fig. 14.41 shows, even at large rapidities the saddle point method 
does not approximate the accurate result reasonably. 




Figure 4.3: The kernel for conformal spin is shown for a s = .1. The differ- 
ent lines correspond to the four resummation schemes, plotted 
in the same style as in Fig. 14.21 



4.1.3 Phenomenology of Mueller-Navelet jets 

Ten years ago the D0 collaboration at the Tevatron measured the azimuthal 
decorrelation between Mueller-Navelet jets |34| . At that time merely the 
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Figure 4.4: The rapidity dependence of the ratio between the saddle point 
approximation and the accurate numeric evaluation for the co- 
efficient Co is shown for a s = .1. The solid line displays the LO 
calculation. The resummation schemes from top to bottom are: 
scheme 2, scheme 1, scheme 3, scheme 4. The lines for scheme 
3 and 4 lie on top of each other. 



LO BFKL equation was available which failed to describe the data by esti- 
mating too much decorrelation, while an exact fixed NLO (af ) Monte Carlo 
calculation by the program JETRAD [103] underestimated the decorrela- 
tion. In contrast, the Monte Carlo program HERWIG |104j was in perfect 
agreement with the data. In Fig. 14.51 we show the data on < cos <p > and 
< cos 2<p > and compare them with our resummed NLO prediction devel- 
oped in the previous section using Eq. (|4.16p . which expresses the angular 
mean values in terms of the coefficients C n . For comparison we give the LO, 
and not-resummed-NLO BFKL results as well. 

It can be seen that the NLO corrections to the BFKL calculation change 
the results significantly leading to a slower decorrelation with increasing ra- 
pidity. As stated in the previous section, the choice of the rapidity variable 
Y = In — ^— turns the convolution with the effective parton distributions into 

PlP2 ^ 

a simple global factor which cancels whenever we study ratios of cross sec- 
tions or coefficients C n . For such observables, the hadronic level calculation 
does therefore not differ from the partonic one. Thereby, uncertainties of 
parton distribution functions do not spoil our calculation. 

It is worth remarking that of all eigenvalues associated to the NLO BFKL 
kernel the one corresponding to the n = conformal spin is the one with 
the poorest convergence in the transition from LO. This becomes evident if 
we consider ratios of coefficients 

C = <«*(m0)> 
C n < cos(n<p) > 

as we did for the case m = 2, n = 1 in the third plot of Fig. 14.51 This 
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Figure 4.5: < cos 4>>,< cos 2</> >, and ^. C C ° S S 2 ^ > at a collider with a cen- 
ter of mass energy ■ v /s=1.8TeV at leading (solid) and next to 
leading order (dashed). The results of the resummation follow- 
ing scheme 1 (dotted) and scheme 3 (dash-dotted) are given as 
well. Tevatron data points are taken from Ref. [105J. 
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observable reveals that the not-resummed-NLO calculation just accidentally 
fits the data for < cos 4> > . 

To further motivate the necessity of a resummed kernel, we now discuss 
the sensitivity against a change of the renormalization scheme. In Fig. I4.6l we 
compare the results for MS renormalization scheme with those for the gluon- 
bremsstrahlung (GB) renormalization scheme |1U01 HOfj . The relevance of 
this more physical renormalization scheme has been discussed in a variety of 
works [601 |6"T| [TOo] 11071 1108j dealing with soft gluon emission. Instead of the 
LO calculation being more dependent on the renormalization scheme than 
the NLO one, nearly the opposite observation has to be made. In contrast, 
the resummed calculations are nearly independent of the renormalization 
scheme, with scheme 3 being the most robust. This underlines the necessity 
to resum such formally sub-leading terms although a quick glance at just the 
observable < cos0 > might suggest that the pure NLO kernel is sufficient. 





Y 




Y 



Figure 4.6: < cos (f> > at a pp collider with a center of mass energy 
y / i=1.8TeV. Results for MS-renormalization are shown as a solid 
line, GB-renormalization is given as a dashed line. The differ- 
ent plots show in the first row LO and NLO, while the second 
row contains the resummation schemes 1 and 3. Tevatron data 
points are taken from Ref. [105J. 



In addition to the deduced observable C m /C n defined in Eq. (|4.33j) . we 
study directly the differential angular distribution as given in Eq. (|4.14p . 
The D0 collaboration published their measurement of the normalized an- 
gular distribution for different rapidities |34j . In Fig. 14.71 we compare this 
measurement with the predictions obtained in our approach from LO, NLO, 
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and resummed BFKL kernel. This comparison throws more light on the 
question whether a resummed NLO kernel is needed. While the weights 
< cos(n0) > seemed to favor the pure NLO calculation, the shape of the 
differential distribution puts this impression into the right perspective. Al- 
though neither the resummed nor the pure NLO kernel matches the data, 
the resummed kernel agrees in shape and is closer to the data than the LO 
estimate. Nevertheless, it still overestimates the decorrelation, and a pure 
calculation of x 2 / n -d-f- still favors the NLO calculation. 

A very strong hint, that the resummation is necessary, is given by the 
instability of the pure NLO kernel against a change of the renormalization 
scheme. Switching to GB-renormalization the angular distribution even be- 
comes negative for large (p. This unphysical behavior clearly is not accept- 
able. 

To comment on the differences between scheme 1 and 3, we have to 
anticipate some details from appendix [B] where we explicate why scheme 3 
is more robust than scheme 1. 

Before we turn to the study of electron-proton collisions in the next 
section, we like to state that we strongly suggest to measure the angular 
correlation of Mueller-Navelet jets at the forthcoming Large Hadron Collider 
(LHC) as well. Besides the possibility of an independent check of the data 
published by the D0 collaboration, measurements for even larger rapidity 
are desirable. This region of forward physics is interesting in many aspects, 
as e.g. to study diffraction or saturation. Since the available rapidity range 
is restricted rather by the geometry of the detector than by the energy 
of the colliding particles, the experimentalists spare no effort to extend the 
capability of the detectors to the very forward region. Therefore, we provide 
our numerical calculations of ratios C m /C n for a broader range of rapidity as 
a prediction for LHC in Fig. 14.81 For the differential cross section itself we 
provide predictions in Fig. 14.91 starting with rapidity Y = 3. Since BFKL 
dynamics require large rapidities, predictions based on the BFKL equation 
for even smaller rapidities are not reasonable. 

Although we soundly justified to ignore the NLO contribution to the jet 
vertex, we studied the possible effects of corrections to it. We investigated 
the impact of the running of the coupling and of that part of the NLO 
Mueller-Navelet jet vertex originating from the splitting functions. These 
contributions can be read off easily from the results of Refs. [251 126] . It 
turns out, that the effect on the overall normalization can become large, as 
it has been shown in Ref. [82J as well, but the change of the ratios that we 
consider is only of the order of a few percent. 

Furthermore, are aware of the fact, that the change from physical ra- 
pidity Ay to an effective rapidity Y is note compensated (See Eqs. (|4.17( 
I4.18P ). We take the resulting uncertainty into account by varying our choice 
of sq = P1P2 by a factor of 2. We varied the renormalization scale fi by a 
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Figure 4.7: ^^jjr at a collider with a center of mass energy y / s=1.8TeV at 
leading (stars) and next to leading order (squares). The results 
of the resummation following scheme 1 (circles) and scheme 3 
(triangles) are given as well. Plots are shown for Y = 1, Y = 3, 
and Y = 5. Tevatron data points with error bars are taken from 
Ref. [33|. 
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factor of two as well, and give the uncertainty arising from these two sources 
by gray bands in Figs. 14.81 and 14.91 




Figure 4.8: Different ratios of coefficients C n obtained by resummation 
scheme 3. The gray band reflects the uncertainty in so and 
in the renormalization scale \i. 
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Figure 4.9: j^^t obtained by resummation scheme 3 for rapidities Y = 
3, 5, 7 from top to bottom. The gray band reflects the uncer- 
tainty in sq and in the renormalization scale fx. 
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4.2 Forward jets at an ep collider 

At an electron proton collision one can study the angular correlation between 
the electron and a forward jet form the proton. A LO study can be found in 
Ref. [109] . Here we apply our NLO formalism to this subject. While the jet 
vertex of the forward jet, in principle, remains the same as in the Mueller- 
Navelet case, we have to deal with a new vertex for the leptonic part, namely, 
for the coupling of the electron to the gluon ladder. Since the electron does 
not interact strongly, the coupling is mediated by an additionally produced 
quark-antiquark pair which couples to the electron via a virtual photon. 
Again we restrict ourselves to LO jet vertices, especially since the challenging 
calculation of the NLO photon impact factor which is needed for the leptonic 
part is still work in progress [HH1 ESI EH ED Mi US] • 

The kinematical situation is depicted in Fig l4.10l where we indicate that, 
for the leptonic vertex, we take the quark-antiquark pair as being inclusive. 
We focus on the outgoing electron with momentum k\ and the gluon with 
momentum q\ which couples to the gluon Green's function. Again, we denote 
the azimuthal angle of k\ as ai and the azimuthal angle of q± as 9\. 




Figure 4.10: Kinematics of the partonic cross section in the ep case 



It is convenient to invoke the common language of deep inelastic scat- 
tering, i.e. denoting the proton momentum as P and the momentum of the 
virtual photon as q 7 , we introduce the virtuality Q = —q^, the Bjorken seal- 
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ing variable xb] = and the inelasticity y = p ^ q q + kl ^ ■ The longitudinal 
momentum fraction of the forward jet will be denoted as xfj- Using the 
relation k 2 = (1 — y)Q 2 and the specific structure of the leptonic vertex, we 
can transform the integration over ki and write 



Hs) =^ J d 2 ki J d 2 k 2 J ^(k 1 |i lcptonic / w l» jcti2 |k 2 )e- y (4.34) 
j d ai j dy j d 2 k 2 j^j d 2 qi J d 2 q 2 J dv J dv' 



2 

n,n'=— oo 



(y,ai\ ^leptonic I Qi ) (qi I v, n) (n , v I f u I v , n ) (n , v | q 2 ) (q 2 1 $j e t,2 1 k 2 ) e^ y , 

(4.35) 



where the rapidity is defined as Y = In xpj / xbj ■ 



4.2.1 The leptonic vertex and a modified jet vertex 

It has been shown in Ref. [109] (taking into account the different conven- 
tions) that the leptonic vertex reads 



(y, ai| Peptonic | qi) 

4a 2 f 1 f 1 <K<k 



2 * u 

l 



7T 2 NcyQ 2z^^j Q Jo e(i-0Q 2 + C(i-C)q! 



x|(j-«i-0- C(i - + 2£(i - 0C(i - C)J y 2 

+ (i - 2^(1 - o - 2C(i - o + - e)c(i - o) (i - v) 

- 4£(1 - 0C(1 - 0(1 - 2/) cos {2(9, - ai)) j 
[ dQ 2 \2af\c L ly,Q 2 ) + 2af\c L ly,Q 2 )cos{2(0 1 -a 1 ))} , (4.36) 



where a denotes the electromagnetic fine structure constant and e 2 is 
the sum over the electric charges of the quark of the quark-antiquark pairs 
which are produced. As a next step we like to compute the projection onto 
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the basis | v, n) . We have 

J d 2 qi (y, ai | $ieptonic I qi ) (qi I v, n) 
= JdQ*J ^! [2aS 0) (q?, y ,Q 2 )+2aS 2) (q?, y ,Q 2 )cos {2(9 1 - a,)) 

= I dQ 2 [2^S 0) (z,,y,Q 2 ) +4 2) (^,y,Q 2 ) (<5 n ,_ 2 e- 2 ' ai + 5 n , 2 e 2iai )" , 

(4.37) 

with 

4 0/2) (,, y ,Q 2 ) = i=|dq? aS°/ 2) (q?,y,Q 2 ) (q?f 4 ■ (4-38) 
To calculate these coefficients, we need integrals of the following type: 



/' 

Jo 



xsri + t f -ii/,i + t f - iiA , (4.39) 

with £? the Euler beta function. Using this formula with £ {0, 1} we 
obtain 



yl\ c H \bu(v z + 1) cosh(7r^j 



V + (12i/ + ll)(l-y) 1 , (4.40a) 

' yN c ' ' H \bu(v z + 1) cosh(7r^) 

x (-(4^ 2 + l)(l-y)). (4.40b) 

Unfortunately, the angular correlation between the electron and a for- 
ward jet has not been measured so far. Nevertheless, to provide theoretical 
predictions ready for a comparison with experiment, we try to implement 
the same kinematical cuts and constraints as the ones that will be used from 
experimental side. Therefore, we consider a vertex for the forward jet, which 
is slightly modified, although we could - in principle - proceed with the jet 
vertex of Eq. (|4.2f) presented in the Mueller-Navelet scenario. The ZEUS 
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collaboration (see e.g. Ref. [80]) and the HI collaboration (see e.g. Ref. |81| ) 
impose also an upper cut for the transverse momentum of the forward jet 
in terms of the photon virtuality Q 2 to ensure that both ends of the gluon 
ladder have a similar scale: 



ZEUS : - < -f < 2 (4.41a) 

HI : \ < M < 5 . (4.41b) 



This requirement suppresses DGLAP evolution without effecting the BFKL 
dynamics. The necessary modification of Eq. (|4.10p is straight forward and 
yields for the ZEUS condition 



c 2 (z/) 



V2 \ + iv' 



(4.42) 



For the HI condition the 1/4 has to be replaced be a 1/10. For simplicity, 
we keep the ZEUS condition. 



4.2.2 Phenomenology of forward jets at an ep collider 



Inserting the leptonic vertex of Eq. (I4.39|) in the leptonic cross section, 
started in Eq. (|4.35p . we get by performing the same transformations as 
in the Mueller-Navelet case 



da 



dy dQ z 



dv I dv' 



+ 4 2) M,Q 2 ) <2,i/|/|z/,2)c 2 (*/)cos 



4 0) {v,y,Q 2 ) (0,H/V,0)c 2 (i/) 



(4.43) 



where we introduced the azimuthal angle 4> = a 2 — ot\ between the elec- 
tron and the forward jet. Furthermore, we made use of the fact that 
(n,v\f\v',n) = (-n,u\f\u',-n). 

We handle the v derivative in f u as we did in Eqs. (|4.23l - 14.24)) . The 
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required terms are 



.d. Af (u,y,Q 2 ) 

z 7r ln r~\ 



21nQ 



+ 



31n(2) 



7T- 



2v 2 + | 5 - 4cos(i/ In 4) 
cosh(2vr^) - 3 4 sin(f In 4) ln(2) 



+ 



8v 



sinh(27r^) 5 — 4cos(f In 4) 
y 2 + 6(l-y) 



+ 



9y 2 + 22(1 - y) + 4(y 2 + 6(1 - y))^ 2 
3^ 2 + l ^ 



v(v 2 + 1) z^2 + i 



(4.44a) 



*7T ln 



21nQ 



+ 



31n(2) 



7T- 



+ 



1v 2 + \ 5-4cos(^ln4) 

cos1i(2tt^) - 3 4 sin(i> In 4) ln(2) 

sinh(27rz/) 5 — 4cos(f In 4) 

3u 2 + 1 3u 



v{y 2 + 1) ^ + l 



(4.44b) 



Hence we can rewrite Eq. (|4.43h as 



da(Y) 



ir 2 a 2 



5 (0) (y, Q 2 , Y) + 5 (2) (y, Q 2 , F) cos 2 J , (4.45) 



# dy dQ 2 2 
where the coefficients are given at LO as 

Bffi (y,Q 2 ,F) = | (^,y,Q 2 ) c 2 (,)e 

and at NLO as 



a s (Q^Y(x (\n\,u)+a s (Q^( X i(\nM-^- cX0 



Ya 3 xo(\n\,v) 



(4.46) 



x e 



21nQ 2 +i 7 fln V 7 



c 2<» 



(4.47) 



For the coefficients with resummed kernel we can write accordingly 

M^) = ^ W (*^ 2 ) c 2 (,)e^- rcsum (l"l^). (4.48) 



When we introduced the new modified jet vertex in Eq. (|4.42p . we were 
considerate of the kinematical cuts applied in experiment for the forward jet. 
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In addition, we use the following cuts motivated by experimental receivables 
[110J concerning the leptonic part: 



20GeV 2 < Q 2 < lOOGeV 2 , .05 < y < .7 . 

5 • 10~ 3 < x Bj < 4 • 10~ 4 . 

Finally, the cross section at the hadronic level reads 

da 



(4.49) 



dY 



: C (Y) + C 2 (Y) cos 2(f>, (4.50) 



Y 



with 

Cn{Y) = ^- f dx F3 dQ 2 dy f eS (x F j,Q 2 )B^(y,Q 2 ,Y)5 (x F3 -^e 

2 ./cuts V V s 

(4.51) 

where we performed the convolution with the effective parton distribution 
from Eq. (|4.5p . By the index at the integral sign, we indicate that the cuts 
of Eqs. (|4.49p are applied. The integration over the longitudinal momentum 
fraction xfj of the forward jet involves a delta function fixing the rapidity 
Y = lnxFj/xBj. Let us note that an additional experimental upper cut on 
xfj would change the single coefficients C n , while the change of their ratios 
is negligible. 

Since the structure of the electron vertex singles out the components 
with conformal spin and 2, the number of azimuthal observables is lim- 
ited when compared to the Mueller-Navelet scenario. Nevertheless, we can 
calculate the rapidity dependence of < cos 20 >= C2/C0. The result is 
shown in Fig. 14.111 and again the NLO calculation predicts a slower decor- 
relation compared to the LO calculation. Since the convolution with the 
effective parton density does not give just a global pre- factor, we have to 
numerically evaluate the kernel several times. For the resummed kernel this 
requires to solve a transcendental equation which turns out to be rather 
time consuming. In the Mueller-Navelet case we recalculated all our results 
with an approximation to the resummed kernel developed in Ref. [60J. The 
agreement is very good, while the time to perform the numerical evaluation 
is drastically reduced. In addition, it approximates the kernel with impact 
factor contribution in the DIS case very well. On that account we present 
the results for this approximation as being representative for the resummed 
kernels in Fig. 14.111 

Because the evolution in rapidity is driven by the kernel, the Y depen- 
dence of the different results is very similar to the Mueller-Navelet case. 
This observation gives, a posteriori, an additional legitimation to consider 
the jet vertices at LO when studying the rapidity dependence. The addi- 
tional inclusive quark-antiquark pair - produced to couple the electron to 
the gluon evolution - yields already in the case of no gluon emission some 
angular decorrelation between the forward jet and the electron. Following 
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Figure 4.11: < cos 20 > at the ep collider HERA at leading (solid), next to 
leading order (dashed) , and for resummed kernel (dash-dotted) . 



the treatment of the Mueller-Navelet case, we present as our final prediction 
the calculation based on the resummed kernel with an error band reflecting 
the uncertainty in so and in \i in Fig. 14.121 For this purpose we varied sq 
and [i independently by a factor of 1/2 and 2 respectively. 
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Figure 4.12: < cos 2<p > at the ep collider HERA for resummed kernel (dash- 
dotted). The gray band reflects the uncertainty in sq and in 
the renormalization scale \i. 



Chapter 5 

Summary and Outlook 



In this thesis we have investigated phenomenological implications of pertur- 
bative QCD calculations on jet observables. We did so in the high energy 
regime, described by the BFKL equation, where our special interest lay in 
effects of NLO BFKL accuracy on exclusive measurements. 

The first part of this work dealt with the production of a jet at central 
rapidity Based on a careful study of the different contributions to the NLO 
BFKL equation and the different energy scales involved, we have derived a 
NLO jet vertex in /c^-factorization. Therefore, we 'opened' the NLO BFKL 
kernel and introduced a jet definition in a consistent way. We have shown 
the infrared finiteness of this jet vertex and its dependence on the scale s\, 
which separates MRK from QMRK. Furthermore, we have examined the 
connections between the energy scale sq, the Reggeon scales SR-ij, and the 
final form of the kernel. 

As a central result of this part, we have calculated the jet production 
vertex (|3.42|) . and we have explicitly given the necessary subtraction, both 
at the matrix element level (|3.40p as well as integrated over the correspond- 
ing phase space (I3.72p , Additionally, we have given an alternative form of 
this subtraction term (|A,41 IA,12j) . It is not as compact as our original one, 
but it combines the divergent parts in a way, that it is very convenient if one 
wants to switch from MS renormalization to gluon bremsstrahlung renor- 
malization. The difference between the two subtraction terms, of course, 
concerns only the finite part of them. Our calculations also suggest that the 
natural scale for the running of the coupling at the jet vertex is the square 
of the transverse momentum of the jet (|3.43p . 

Our vertex can be used in the context of 7* 7* -scattering (|3.4p or of 
hadron-hadron scattering (|3.59|) . The original derivation of the BFKL equa- 
tion, at LO and NLO as well, relies on the colliding objects to have a hard 
and similar scale. The scattering of virtual photons is a natural implemen- 
tation of this requirement. Hence, for this process the embedding of our 
jet vertex is straight forward. We have shown that for hadron-hadron scat- 
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tering special care has to be taken. The asymmetry between the soft scale 
of the hadrons and the hard scale of the jet induces additional evolution 
in kx- This fact is reflected by a modified kernel (13.531) which governs the 
evolution of the unintegrated gluon density, and a modification of the jet 
vertex itself (|3.54p . The modification to the evolution kernel can be consid- 
ered as the first term of a resummation of higher order terms, induced by 
the scale change. Hence, for the first time, we have given an example of a 
complete process described in NLO /cr-factorization with NLO integrated 
gluon densities. Nevertheless, we have to state that our expression is only 
valid for small x. 

In our analysis we have been careful to retain the dependence upon the 
energy scale sa which appears at NLO accuracy and separates multi-Regge 
kinematics from quasi-multi-Regge kinematics. In the NLO calculation of 
the total cross section, one may be tempted to take the limit s\ —* oo, thus 
disregarding the 1/sa corrections to the NLO BFKL kernel. However, when 
discussing inclusive (multi-) jet production one has to remember that s\ has 
a concrete physical meaning: it denotes the lower cutoff of rapidity gaps and 
thus directly enters the rapidity distribution of multi-jet final states. In a 
self-consistent description then also the evolution of the unintegrated gluon 
density has to depend upon this scale. 

In fact, to study multi-jet production is a promising application of our 
jet vertex. The vertex, as we provide it, is already well suited for the im- 
plementation into a Monte-Carlo program. There, momentum conservation 
can be taken into account, and together with the running of the coupling 
and a physical choice of s\ deeper insight into QCD multi-jet production is 
expected. Especially, the large available center of mass energy at the forth- 
coming Large Hadron Collider (LHC) calls for an appropriate description 
of multi-jet events. The LHC will only tap its full discovery potential if 
a well understanding of the background is guaranteed. Furthermore, the 
dependence on s\ is an inherently interesting topic. Therefore, a numerical 
study of this feature will be advisable. For this purpose, the calculation 
of the inclusive jet cross section, as presented in our work, might act as a 
theoretical laboratory environment. 

Although the NLO unintegrated gluon density might appear as some 
kind of byproduct of our work, it is a significant result. Further applications 
of our NLO /c^-formalism include W and Z as well as heavy flavor production 
in the small- x region. Compared to the results presented in this paper, these 
applications require the calculation of further production vertices; however, 
for the treatment of the different scales and of the unintegrated gluon density 
all basic ingredients have been collected in this thesis. 

In the second part of this thesis, we studied the angular decorrelation 
between Mueller-Navelet jets, and between electron and forward jet in deep 
inelastic scattering (DIS). We have derived an analytical LO master formula 
(14.141 14~4"5|) . with only one integration left to be performed numerically. Fur- 
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thermore, we have obtained the corresponding expression including the NLO 
kernel in its angular dependent form (|4.28|) . In face of the scale dependent 
part, we have shown, that the kernel still can be exponentiated from the very 
beginning. In our approach we circumvent to be content with an expansion 
of the Green's function. As a consequence, a contribution from the impact 
factors entered the energy dependent exponent. This indicates that the full 
separation of energy dependence and impact factors does only hold at LO 
and is violated due to the running of the coupling, spoiling the conformal 
invariance. 

The NLO corrections to the kernel are known to be negative and very 
large. It has been shown that it even leads to instabilities which can be 
cured by a resummation of terms beyond NLO. Therefore, in section 12.31 
we extended different resummation procedures to arbitrary conformal spin. 
We then derived the equations describing angular decorrelation also for the 
resummed kernel (|4.29|) . We explicitly demonstrated the need for this re- 
summation by studying the dependence on renormalization schemes and 
scales. 

We have compared our results for Mueller-Navelet jets with data ob- 
tained at the Tevatron. The comparison is given in Figs. 14.51 14.71 We have 
reproduced the LO result that overestimates the decorrelation. In contrast, 
our calculation using the resummed kernels improves the theoretical predic- 
tion although it still overestimates the decorrelation. Our study of different 
renormalization schemes has given additional support to the necessity of re- 
suming the kernel. We strongly suggest to study the different observables, 
connected to angular decorrelation of Mueller-Navelet jets, at the LHC as 
well. For that purpose, we give predictions in Fig. 14.81 where we also have 
elaborated the theoretical uncertainty by the renormalization scale \i and 
the energy scale so- 

Finally we have extended our theoretical work to the scenario of DIS 
where one can study the azimuthal correlation of the electron and a forward 
jet. Unfortunately, there are no corresponding measurements. None of the 
forward jet studies at HERA has studied the angular dependence. We hope 
that the predictions we give - incorporating the experimental cuts usually 
applied in studies of forward jets - will motivate the experimentalists to 
investigate this observable. 



Appendix A 

Alternative subtraction term 



In this appendix we present an alternative subtraction term which does not 
make use of the simplifications A^ + — > 2A^ and + A^ — > 2A^ 
which we used in Eqs. (|3.36l I3.38|) . These limits are valid in the kinematic 
regions leading to IR-divergences and hence they do provide the correct e 
poles. However, they also alter the finite terms. Here we want to study 
also this finite part as accurately as possible and hence we do not take these 
limits but use instead the complete sum 

A {1) + A {2) + A (3) + A (4) + A {5) + A (6) + A M rk (A.l) 

as the gluonic subtraction term. 

The full gluonic matrix element written in Eq. (|3.27|) contains spurious 
UV-divergences which are canceled when combined with the MRK contri- 
bution. One fourth of the MRK contribution cancels the UV-divergence of 
A( 4 ) while another fourth cancels that of A^ . The remaining half cancels 
the UV-divergence of two terms present in Eq. (13.270 : 

A _q'q6 (i-x 1 x i\ 

A ^ =^e~ l^kf + — kf ) ' (A - 3) 

which are IR-finite and hence so far not included in the subtraction term. 

By doubling Au\ and A^ in the subtraction term constructed in Eq. (|3.4U[) 
also their spurious UV-divergences are doubled and thus completely can- 
celed by the MRK contribution. But Eq. (|A.ip so far only contains half of 
the spurious UV-divergences of the full matrix element in such a way that 
half of the MRK contribution is not compensated. Therefore a subtraction 
term based on Eq. (jA.lh which is also free from spurious UV-divergences 
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should also include A™ and A($\ and reads 



H8) 



^uonf = A {1) + A (2) + ^(3) + A( 4 ) + A (5) + A (6) + A M RK + %) + A (8) 

= + A (2) + A (3) + (A (5) - A (4) ) + A {6) + + A (7) + A {8) . (A.4) 



If we now define 5, 



(3,6,7,8) 



and 5mrk as we did in Eq. (|3.6ip we get a new 



integrated subtraction term from the previous Eq. (|3.72p by replacing 



Sin + Siv 



1 tiT(1 - e) / A 



A 2 (4vr) e 



.2\ e 



I* 



1 5vr 2 „ . , 

3-T- + 0(£) 



(A.5) 



with 



i (5/// + 5jy ) + 5(3) + 5 (6) + C?M f K + 5( 7 ) + 5 {8 ) . (A.6) 



The results for 5( 3 ) and Sr$\ can be easily obtained from Eqs. (C.43) and 
(C.40) of Ref. [ITT] : 



'(3) 



'(6) 



1 vrr(l-e) ^ A 2 v 



1 



A 2 (4vr) e v A 4 

q^A^-qj,)) j 



2e 2 2e A 4 12 4 q 2 



+ 



In 



q|qjA^ 

A 2 (q a -q b ) 2 \2~\dU~\(<& + <$) 4 



- Li 2 



q| 



+ Li 2 



In 



dz 



q« 



q|q|/ (A(q a -q fe ))' 



A 2 (q a - q 6 ) 5 



(q 6 + zq a ) 2 



+ 0(e) 



1 vrr(l - e) / A 



(47T) 



.2\ e 



1 7T 



e z 6 



(A.7) 
(A.8) 



Due to the UV-singularity of Amrk we regularize the x integration by a 
cutoff 5 to obtain 



5mrk 



1-5 



dx 



d D ~ 2 k 7 



1 



(l-x)J /i 2e (2^- 4 k 2 (A-k 2 ) 2 
1 Trr(l-e) ^A 2 Vr 2 (6) 2]n 5 



(47T) 



r(2e) 



1-5 



(A.9) 



Making use of 2q a ki — q 2 = t + k 2 /x we can decompose Eq. (C.41) of 
Ref. [Ill] into one integration very similar to that of 5mrk and another one 
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which can be transformed to give 5(7). 



>(7) 



1 tiT(1 - e) 
A 2 (4vr) e 



A 2 



1 r 2 (6) 

2 r(2e) "'1-5 



In- 



1 
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1 n 2 7T 2 
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The two parts forming can be obtained from each other by the exchange 
k% <-» A;2 and we only need to double the calculation of one: 



>(8) 
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dx 



x(l — x) 

l ~ S dx 



d D - 2 k, 



X 
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2T(2e) 

When we add up these new contributions, the spurious UV-divergences in- 
deed cancel and we can safely take the 5^0 limit. Furthermore, the 
new subtraction term has the same pole structure and differs only in the 
finite parts when compared to that in Eq. (|3.40p and its integrated form in 
Eq. (|3.72p . To complete the calculation we combine it with the correspond- 
ing unmodified quark part and obtain 



d D ~ 2 k 2 
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Appendix B 

Resummation of sub-leading 
corrections to the BFKL 
kernel 



In this appendix we want to elaborate the different schemes of resummation 
of formally sub-leading corrections to the BFKL kernel, as it is introduced 
in section 12.3.21 and used in chapter 01 

The first section describes different concrete schemes to implement the 
resummation procedure. We then briefly show how a different renormaliza- 
tion scheme effects the resummation. We have shown in chapter H] how the 
running of the coupling effectively leads to an contribution of the impact 
factors to the kernel (see Eqs. (14. 19111311) ). In section fB.31 we explain how 
the resummation changes if this contribution is included. In the last section 
of this appendix we discuss some consequences of the resummation in detail. 

B.l Different schemes of implementations 

B.l.l Scheme 1 

For scheme 1 we choose £^(7) to be 

Dh{l) = { {k -1)1 ^ ~ ^ (7)] - 

From Eq. (|2.58j) we get 

x (°) (m, 7 ) = 2V(1) - V>(7+Y + !)-^(l-7+y + |). (B-2) 

By an cj-expansion of Eq. (|B.2h we obtain 

(o), 1 r ,/ / m\ , / m\~\ 

xi \™>v = -2^0(^,7) \w [1 + +v ^1-7 + yJJ • ( B - 3 ) 
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From Eq. (|2.56j) we calculate d^\{m) = dr^{m) an< ^ 



df[(m) = -*'(l + m), 



7(0) 
1,1 

M 
n,2 

7(0) 



d$M= -^(l + m)-*(l) 
1 

~ 2' 



d$(m) = 

We can now use master formula (|2.59p to write 

(771,7) = X (0) (w-,7) +«s (xi(jn,-y) - Xi 0) ("7,7)) 
2 

+ a s ^2 ( d i,k( m ) ~ d i,k( m 



(B.4a) 
(B.4b) 

(B.4c) 



k=l 



Dk ^ + J + 2)- Dk b + J 



777 UJ 



+ D k [l- 1 + - + -)-D k [l- 1 +- 



m 



. (B.5) 



The sum in Eq. (1B.5|) can be truncated after k = 2 since (74,3(777) = 
B.1.2 Scheme 2 

For scheme 2 we choose a simpler £>/<; namely 



in 



1 



(B.6) 



From Eq. (12.580 we get 

X (0) (t77,7) = Xo("7,7) 



1 



7 +f 1-7+f 
1 

+ 



+ 



, (B.7) 



and 



Xi ("7,7) 



Xo(m,7) 



1 



+ 



1 



(7 + yj- (l-7 + §)7 
Also for this choice we have d^l(m) = <l^l(m) which explicitly read 



(B.8) 



"1,1 



<42- ll' 1 '" + *<•>)• 



#'(l+m) - tf'(l) + 



1 



(1 + 777) 



(B.9a) 

(B.9b) 
(B.9c) 
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B.l. 3 Scheme 3 

Scheme 3 is based on scheme 1. It is then modified to hold d^ k (m) = d\^(m) 
not just for k = 3, but also for k = 1, 2. To do so, we modify Eq. (|B.2p by 
two terms of higher order labeled A(m) and B{m): 



x (°)( m , 7 ) =(l-a 8 A(m)) 



2iP(l)-ij[ 1 + - + - + a s B{m] 



, m u) „, . 

^(1-7 + ^+ 2+ OLsB{m) 



(B.10) 



From w-expansion we get: 



Xi ( m ,V 



i \ r / Tn \ / Tft 

B ( m ) + ^XoO™, 7) j |V> (7 + + ip [1 - 7 + 2" 



A(m)xo{rn,l)- 



(B.ll) 



The coefficients dS ' now depend on these newly introduced variables A and 
B in the following way 



dg[ (m) 
djg(m) 

dg(m) 



^4(m) - + m), 

B(m) + -(^(l + m) - 
1 

2' 



(B.12a) 
(B.12b) 

(B.12c) 



The requirements d^\{m) = d 1; i(m) and d^°](m) = di,i(m) can be satisfied 
by an appropriate choice for A(m) and B(m): 



A(m) =-di,i(m) -*'(l+m) (B.13a) 
S(m) = - di )2 (m) + - (^(1 + m) - * (1)) . (B.13b) 
Due to that extra requirement the master formula simply reads 

X (1) (m,7) = x (0) («i,7) +"s (xi(™,7) ~ Xi 0) ("i,7)) • ( B - 14 ) 

B.l. 4 Scheme 4 

Scheme 4 uses the same functions as scheme 2, but is modified similar to 
scheme 3 to hold d^\{m) = d\ k(m) not just for k = 3, but also for k = 1,2: 

1 1 



X (0) ( m >7) = Xo(m,7) 



1 



7+f 1-7+f 



+ 



+ (1 - a,A(m)) 



1 



7 + § + % + a s B{m) 1-7 + f + f + a s S(m) 



. (B.15) 
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From w-expansion we obtain: 



(°)/ ^ 



(1 + m)A(m) 



7 + f 1-7+f 



5(m) + 2^0(^,7) 



+ 



(7 + f) 2 (1-7 + f) 2 . 



The result of calculating the coefficients d^ ' is 



(B.16) 



d?hm)=-A(m)-~ 



V(l+m) - *'(1) + 



(1 + m,y 



dg(m) 

(45 m 



B(m) + -(¥(l + m)- 
1 
2 



(B.17a) 

(B.17b) 
(B.17c) 



The requirements d^\(m) = dii(m) and d^\{m) = d\ \(rn) can be satisfied 
by an appropriate choice for A{m) and B{m): 



A(m) = — d\ : \(rn) 



*'(l + m) - *'(1) + 



1 



(1 + mf 



B(m) 



di,2(m) + -(*(! + m) - 



(B.18a) 
(B.18b) 



The resulting master formula equals that of scheme 3 given in Eq. (|B,14p . 
but with the different x^ \ m > 7)> Xi (^jT)- 



B.2 Renormalization schemes 

To perform the renormalization in the gluon-bremsstrahlung (GB) scheme 
instead of the MS scheme is straight forward. The redefinition of the Landau 
pole 

2Nc_ <, 

A m - A GB = Ams e "0 (B.19) 

is accompanied by canceling the 5-term in Xi(m, 7) and, as a consequence, 
in o?i ) i(m) as well. 

B.3 Impact factor contribution 

Including the impact factor contribution of Eq. (|4.25h to the kernel in the 
Mueller-Navelet case effectively changes XV 

Xi(m, 7 ) - Xi(m, 7 )-^- 7(1 _ 7) , (B.20) 
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leading to a change in d\^(m): 

I / \ , / x . A) ( 1 — S m ,0 Sm,o\ Cant \ 

d ltl (m) -> dl>1 ( m ) + _^_^_^-_j (B.21a) 

di >2 (m) -> di )2 (m) - wrr^mfi- (B.21b) 

oiv c 

In the DIS case the the additional contribution of Eq. (|4.44j) effectively 
modifies xi an d hence the coefficients di^(m): 

44(0) - ^ l( 0)-A^ + y(| l_-» + j (B.22a) 

rfi, 2 M di, 2 (m) + ^r. (B.22c) 



B.4 Discussion 



In Fig. IB.ll we show the kernel for conformal spin m = to m = 2 in 
dependence on 7. An important common feature is, that the the LO kernel 
has simple poles at 7 = — m/2 and 7 = l + m/2, and the NLO kernel simple, 
double, cubic poles at these points. After the resummation procedure these 
poles are removed. For conformal spin m = 2, we observe simple poles 
at 7 = and 7 = 1 for the NLO and resummed kernel, originating from 
the term 

^ ~ — 77 — \^7ti 2 h^ the NLO kernel (see Ecj. f|2.46p 1 . Similar poles 
arise for all conformal spins m at 7 = and 7 = 1 when the impact factor 
contribution, shown in the previous section, is included. These are not 
resummed as well. The large difference between LO and NLO kernel around 
7 = 1/2 is specific to m = 0. Therefore, also the resummed does not deviate 
visibly from LO and NLO kernel around 7 = 1/2 for m > 0. This indicates 
that a resummation of the poles at 7 = and 7 = 1 for higher conformal 
spin is not imperative. 

In Fig. IB.2l we show the kernel for conformal spin m = 0. The completely 
different shape around v = for the LO and NLO version of the kernel can 
be clearly seen. Since the results of the resummation scheme 2 (4) is hardly 
distinguishable from scheme 1 (3), we focus on schemes 1 and 3 for the 
plots. In deed, the resummation procedure leads to one maximum at v = 
and agrees with the NLO results for larger \u\. A second striking feature of 
the resummed kernel is its stability against a change of the renormalization 
scheme. We have checked that for larger conformal spin LO, NLO and 
resummed kernels behave in the same way. For this reason we do not present 
the according plots. 
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Figure B.l: Comparison of kernels uj{m = 0,7) at a scale /i = 30GeV as a 
function of 7. Shown are the LO expression a s Xo (solid line), 
NLO expression a s Xo + ®^Xi (dashed line), and the resummed 
one for scheme 3 (dashed-dotted line). 
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w (m = 0, 1 + if) w (m = 0, \ + if) 




Figure B.2: Comparison of kernels uj(m = 0, ^ + iv) at a scale \i = 
30GeV as a function of v. The left diagram shows the re- 
sults in MS-renormalization and the right one the results in GB- 
renormalization. Shown are the LO expression (solid line), NLO 
expression (dashed line), and the resummed ones for scheme 1 
(dotted line) and scheme 3 (dashed-dotted line). 
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Appendix B. Resummation 



When we extend the resummation prescription to the exponent of the 
coefficients C n defined in Eq. (|4.28|) . the kernel is augmented by an impact 
factor contribution ~ Xo/(t(1 — 7)) with its implications shown in section 
IB.31 The shapes of the LO, NLO, and the resummed version of the kernel 
remain the same compared to the pure form without impact factor contri- 
bution (see Fig. IB.3[) , but scheme 1 and 2 are more sensible to the depth of 
the NLO-dip around v = 0, because in these schemes terms such as a"/7 n+1 
have not been resummed [54J. For scheme 1 and to this can even change 
the maximum at v = back to a local minimum. Since schemes 3 and 4 are 
more stable, we consider them to be the best candidates for a meaningful 
resummed kernel. 
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w (m = 0, 1 + if) w (m = 0, \ + if) 




Figure B.3: Comparison of kernels io(m, + w) including the impact factor 
contribution at a scale /i = 30GeV as a function of v. The 
left column shows the results in MS-renormalization and the 
right one the results in GB-renormalization. Shown are the 
NLO expression (dashed), and the resummed ones for scheme 1 
(dotted) and scheme 3 (dashed-dotted). 
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